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Ultrastructural Studies of Oogenesis in the Fern 
Diplopterygium glaucum 


Xi-Linc Dat’, Jia-QuaN Bao’, QuAN-X1 Wanc!?, Yu-CHEN Cao, AND JiAN-Guo Cao!* 
1. College of Life Sciences, Shanghai Normal University, Shanghai200234, China 


2. Shanghai Key Laboratory of Plant Functional Genomics and Resource, Shanghai Chenshan 
Botanical Garden, Shanghai 201602, China 


Asstract.—The present investigation examines the oogenesis of the gleichenioid fern, 
Diplopterygium glaucum (Thunb. ex Houtt.) Nakai using transmission electron microscopy. The 
results demonstrate that the oogenesis of D. glaucum is similar to the core leptosporangiate ferns. 
The egg cell of D. glaucum was progressively isolated from the adjacent cells by forming a 
separation cavity and an egg envelope during maturation. However, during this process, a pore 
region consistently connected the egg and the ventral canal cell. Ultimately, a fertilization pore 
formed at the pore region when the egg matured. The nucleus produced conspicuous evaginations 
during later stages of egg development. The cytological features during oogenesis in D. glaucum 
more Closely resembled the core leptosporangiate ferns rather than the basal leptosporangiate ferns 
(Lygodium and Osmunda), which possess no egg envelope or fertilization pore. The results provide 
new characteristics for interpreting the evolutionary history of the gleichenioid ferns. 


Key Worps.—gleichenioid fern; egg envelope; fertilization pore 


Observations of oogenesis among fern species representing the core 
leptosporangiate ferns including Pteridium aquilinum (Bell and Miihlethaler 
1962a, b; Bell and Duckett, 1976; Cao, Dai, and Wang, 2012b), Histiopteris 
incisa (Bell, 1980), Athyrium filix-femina (Fasciati et al., 1994), Dryopteris 
crassirhizoma (Bao et al., 2005), Ceratopteris richardii (Lopez-Smith and 
Renzaglia, 2008), Ceratopteris thalictroides (Cao, Yang, and Wang, 2009; Cao, 
Wang, and Bao, 2010a), Adiantum flabellulatum (Cao et al., 2010b), Plagiogyria 
euphlebia (Cao, Dai, and Wang, 2011), Coniogramme emeiensis (Wang et al. 
2012a), Cibotium barometz (Wang et al. 2012b), Anisocampium sheareri (Yang 
et al. 2013), and Phymatosorus hainanensis (Zou et al. 2014) demonstrate that 
all mature eggs of these species possess an egg envelope consisting of a 
reticular or multilayered structure. In addition, all of these species possess a 
fertilization pore, a structure that limits penetration of the sperm (Cao, Yang, 
and Wang, 2009; Cao, Wang, and Bao, 2010a; Cao et al., 2010b; Cao, Dai, and 
Wang, 2011; Cao, Dai, and Wang, 2012b; Wang et al., 2012a, b; Yang et al., 
2013, Zou et al., 2014). The structure and formation of the fertilization pore in 
the ferns had been described in Ceratopteris and Pteridium (Cao, Wang, and 
Bao, 2010a; Cao, Dai, and Wang, 2012b). 

The basal groups of the leptosporangiate ferns include osmundaceous ferns, 
filmy ferns, gleichenioid ferns, and schizaeoid ferns (PPG I, 2016; Pryer et al., 
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2004; Smith et al., 2006). Among these, Osmunda, possesses no egg envelope 
or fertilization pore in the mature egg (Cao, Dai, and Wang, 2012a). Lygodium 
japonicum, representing the schizaeoid ferns, only forms a pore region 
between the egg and the ventral canal cell during oogenesis, with no egg 
envelope or fertilization pore (Cao et al., 2017). 

Gleichenioid ferns lie between osmundaceous ferns and schizaeoid ferns in 
the fern phylogeny (PPG I, 2016; Smith et al., 2006). However, the cytological 
features of the oogenesis of gleichenioid ferns have not been investigated. 
Whether the cytological features of the oogenesis of gleichenioid ferns reflects 
their phylogenetic status remains obscure. In the present investigation, 
oogenesis of the fern Diplopterygium glaucum is observed in detail, and the 
results provide some new evidence for understanding the phylogenetic 
position of the gleichenioid ferns. 


MATERIALS AND METHODS 


Spores of Diplopterygium glaucum (Thunb. ex Houtt.) Nakai (before revision 
called: Hicriopteris glauca (Thunb.) Ching) were collected from plants in 
Ningbo city, Zhejiang Province, China. The spores were surface sterilized with 
5% sodium hypochlorite solution for 3 min. After rinsing three times with 
sterile distilled water, the spores were sown on a modified Knop’s solution 
(0.8g Ca(NO3)2°4H2O; 0.2g KH2PO4; 0.2gKNO3; 0.2g MgSO4°7H20, dissolved 
in 1 liter distilled water), solidified with 1.5% agar in culture dishes. These 
dishes were placed in an artificial climate chamber under conditions of 25°C in 
the light (18 h) and 20°C in the dark (6 h). After 11 to 12 weeks, archegonia had 
developed on the lower surface of the gametophytes just behind the growing 
apex. 

Gametophytes bearing various archegonial stages were placed in 3% 
glutaraldehyde in 0.1 mol/L phosphate buffer at room temperature for 6-12 
h. The specimens were subsequently washed three times with the same buffer, 
postfixed in 2% aqueous osmium tetroxide for 2 h, rinsed three times in buffer 
and embedded in Spurr’s resin (SPI-Chem, USA) via a graded acetone series. 
Specimens were thick sectioned for the presence of the archegonia and thin- 
sectioned with a diamond knife on an ultramicrotome (Leica EM UC7, 
Germany). The thin sections were stained with uranyl acetate and lead citrate. 
All specimens were observed with a transmission electron microscope (Tecnai 
G2 Spirit BioTWIN, FEI Company). 


RESULTS 


Young Egg 

Archegonia of Diplopterygium glaucum are produced on the lower surface just 
behind the growing point of the gametophyte. At the time of egg formation, the 
archegonium contains an axial row of three cells, i.e., an egg cell, a ventral 
canal cell (VCC), and a single elongated neck canal cell (NCC) (Fig. 1a). The 
inner three cells are closely appressed to the archegonial jacket cells (Fig. 1a). 


DAI ET AL.;: OOGENESIS IN THE FERN DIPLOPTERYGIUM GLAUCUM 3 


BS cag: 
1 umpes 
am < 


Fic. 1. YOUNG EGG STAGE OF DIPLOPTERYGIUM GLAUCUM. (a) A young archegonium, containing a newly 
formed egg (E), a ventral canal cell (VCC) and a neck canal cell (NCC) with two nuclei. (b) Part of 
the egg and the VCC showing plasmodesmata (pd) between the two cells, abundant Golgi bodies 
(G), vesicles (Ve) in the upper part of the egg cytoplasm. (c) Lower part of the egg showing thick 
wall between the egg and the jacket cell, abundant cytoplasmic organelles. (d) Lateral part of the 
VCC showing abundant Golgi bodies. (e) Lower part of the VCC, showing abundant Golgi bodies 
and other organelles. (f) Upper part of the VCC, showing abundant endoplasmic reticula (ER). M, 
mitochondrion; N, nucleus; P, plastid. 
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There are well developed plasmodesmata between the egg and the VCC (Fig. 
1b), and between the VCC and NCC (Fig. 1e), but no plasmodesmata were 
observed between the inner cells and the jacket cells (Fig. 1c). The nucleus of 
the egg is ellipsoid in shape and usually contains an irregular nucleolus (Fig. 
1a). Vesicles, with a diameter of 0.5-11m, are frequent in the egg cytoplasm 
(Fig. 1a-c). Plastids contain fewer lamellae and no starch grains (Fig. 1c). 
Mitochondria are distributed throughout the cytoplasm of the egg. In the VCC 
and NCC, there are more Golgi bodies in the cytoplasm (Fig.1 d-e). 
Endoplasmic reticula are frequent in the periphery of the NCC cytoplasm 


(Fig.1 f). 


Maturing Egg 
FORMATION OF THE SEPARATION CAVITY AND Pore REGION.—Soon after the egg is 


formed, a separation cavity forms between the egg and the ventral canal cell 
(Fig. 2 a, b). However, a connecting region (pore region) still links the egg and 
the VCC in the central region (Fig. 2a). There are well developed 
plasmodesmata in the pore region (Fig. 2c). An obvious wall, lying closely to 
the VCC, separates the egg and the VCC (Fig. 2a-c). Some amorphous materials 
can be seen in the separation cavity (Fig. 2b). At this stage, the Golgi bodies in 
the egg cytoplasm increase in number especially around the periphery of the 
egg (Fig. 2d). Some of the mitochondria beneath the pore region become 
elongated (Fig. 2c). Plastids are similar to the previous stage. There are many 
vesicles in the cytoplasm of the egg (Fig. 2a, b). Sometimes, vesicles fuse with 
the separation cavity (Fig. 2b). In the canal cells, ER and Golgi bodies are more 
abundant than before (Fig. 2e, f). Plasmodesmata continually connect the VCC 
and NCC during the development of the egg (Fig. 2e). 


THE FORMATION OF THE EGG ENVELOPE AND THE FERTILIZATION Pore.—At the later 
stage of egg development, the VCC separates from the egg in the pore region 
(Fig. 3a). An egg envelope is formed in the upper surface of the egg (Fig. 3a). 
Higher magnification observation shows that the egg envelope is formed by 
opaque materials deposited on the dark stained plasmalemma (Fig. 3b). 
However, no egg envelop is formed in the pore region (Fig. 3b). The only 
membrane covering the fertilization pore is the plasmalemma (Fig. 3b). The 
diameter of the fertilization pore reaches about 3.5 pm (Fig. 3b). Concentric 
membranous vesicles are seen beside the pore region (Fig. 3a, b, arrows). At 
this stage, the nucleus becomes irregular (Fig. 3b). Abundant vesicles and 
Golgi bodies are mainly located in the central part of the egg cytoplasm above 
the nucleus (Fig. 3a, b). 


Mature Egg 

When the egg matures, a layer of egg envelope covers the surface of the egg 
apart from the fertilization pore (Fig. 3c). The egg envelope in the upper 
surface of the egg is especially prominent. The thickness of the upper egg 
envelope decreases from the center to the periphery (Fig. 3c, d). In the lower 
part of the egg, the egg envelope is much thinner compared to the upper layer 
of the egg envelope (Fig. 3c). There are electron-light spaces inside the egg 
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Fic. 2. MATURING EGG STAGE OF DIPLOPTERYGIUM GLAUCUM. (a) A separation cavity (SC) forms above the 
egg. A pore region (PR) connects the egg (E) and the VCC. (b) Higher magnification of Fig. 2a, a cell 
wall (W) lies close to the VCC. Vesicles of the egg present in separation cavity (SC). (ce) 
Magnification of the pore region (PR) showing plasmodesmata (Pd) in the pore region, 
Mitochondria (M) are abundant beneath the pore region. (d) Lateral part of the egg showing 
Golgi bodies. (e) Part of the VCC and NCC showting abundant endoplasmic reticulum (ER). (f) Golgi 
bodies (G) and ER become more abundant in NCC. G, Golgi bodies; M, mitochondrion; N, nucleus; 
S, starch. 
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Fic. 3. MATURING AND MATURE EGG STAGE OF DIPLOPTERYGIUM GLAUCUM. (a) A fertilization pore (FP) and 
egg envelope (EE) form at the late stage of egg (E) development. (b) Magnification of the fertilization 
pore (FP), abundant vesicles (Ve) form in the upper part of the egg; c. A mature egg with a 
fertilization pore (FP). (d) Magnification showing vacuolation of the egg cytoplasm. Nuclear 
evaginations (arrow) form. (e) The nucleus of the egg becomes irregular in outline with nuclear 
evaginations (Nev). (f) Higher magnification showing the egg envelope and nuclear evaginations. 
M, mitochondrion; N, nucleus; VCC, ventral canal cell; Ve, vesicles. 
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envelope, resulting in a network structure in section (Fig. 3d). Higher 
magnification observation shows that the egg envelope is associated closely 
with the dark stained plasmalemma (Fig. 3f). 

The nucleus of the egg remains highly irregular in outline and produces 
numerous sac-like evaginations near the main body of the nucleus (Fig. 3d, 
arrow; Fic. 3e, f). Some evaginations remain connected to the main body of the 
nucleus via narrow isthmuses (Fig. 3f, arrow). Sac-like evaginations, which, in 
this plane of section, are not connected to the nucleus, are frequently seen in 
the cytoplasm of the egg (Fig. 3f, asterisks). The matrix of these evaginations 
resembles that of the main body of the nucleus (Fig. 3f). However, the double 
membranes of the evaginations cannot be identified clearly (Fig. 3f). It is 
notable that the cytoplasm, surrounded by the cup-shaped nucleus in the 
upper part of the egg, is highly vacuolated (Fig. 3d). At this stage, plastids 
contain no starch grains or lamellae. Plastids can only be recognized by the 
presence of the plastoglobuli (Fig. 3f). The mitochondria possess well- 
developed cristae. At the last stage of egg maturation, the canal cells 
degenerate completely. Almost all of the cytoplasm of the canal cells has 
decomposed into amorphous mucilaginous materials (Fig. 3c). 


DISCUSSION 


The Egg Envelope and Fertilization Pore 

Since the discovery of the fertilization pore in Ceratopteris thalictroides (Cao, 
Yang, and Wang, 2009), we have been dedicated to expanding the number of 
species for which egg development has been ultrastructurally studied to 
examine whether the fertilization pore exists in ferns in different evolutionary 
lineages. So far, it has been confirmed that all core-leptosporangiate ferns 
examined possess the egg envelope and fertilization pore in their mature eggs 
(Cao, Yang, and Wang, 2009; Cao, Wang, and Bao, 2010a; Cao et al., 2010b, 
Cao, Dai, and Wang, 2011, 2012b; Wang et al. 2012a, b; Yang et al. 2013; Zou et 
al. 2014). Among the basal leptosporangiate ferns, Osmunda japonica and 
Lygodium japonicum have been investigated. Neither of the two species 
formed an egg envelope or fertilization pore (Cao, Dai, and Wang, 2012a; Cao et 
al., 2017). The present investigation shows that D. glaucum, a gleichenioid 
fern, which current phylogenetic analyses (PPG I, 2016; Shen et al., 2017; 
Smith et al., 2006) demonstrate is a more ancient lineage than L. japonicum, 
has formed a typical egg envelope and fertilization pore. The diameter of the 
fertilization pore of D. glaucum (about 3.5 um) is obviously larger than those of 
core-leptosporangiate ferns, whose fertilization pores are usually about 2.5-3.0 
uum (Cao, Yang, and Wang, 2009; Cao et al., 2010b, Cao, Dai, and Wang, 2011, 
2012b; Wang et al. 2012a,b; Yang et al. 2013; Zou et al. 2013). The larger 
fertilization pore may be an ancestral feature in fern oogenesis. 

The present investigation shows that the egg envelope of D. glaucum is 
reticulate in structure, which is the same as that found in the Cyatheales and 
most members of the Polypodiales, including Pteridium (Bell and Miihlethaler, 
1962a, b; Cao, Dai, and Wang, 2012b), Athyrium filix—femina (Fasciati et al, 
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1994), Dryopteris crassirhizoma (Bao et al, 2005), Plagiogyria euphlebia (Cao, 
Dai, and Wang, 2011), Anisocampium sheareri (Yang et al., 2013), and 
Phymatosorus hainanensis (Zou et al., 2014), but differs from the egg 
envelopes of the family of Pteridaceae sensu lato (Smith et al., 2006), 
including Ceratopteris (Cao, Yang, and Wang, 2009; Cao, Yang, and Bao, 
2010a; Lopez-Smith and Renzaglia, 2008), Adiantum flabellulatum (Cao et al., 
2010b) and Coniogramme emeiensis (Wang et al., 2012a). Members of the 
Pteridaceae sensu lato usually possesses multilayered egg envelopes, which 
are considered to be formed by attachment of the endoplasmic reticulum (Cao, 
Yang, and Wang, 2008; Cao, Wang, and Bao, 2010a; Cao et al., 2010b). The 
formation of the egg envelope and fertilization pore appears to prevent 
polyspermy (Cao et al. 2010c). However, the two types of the egg envelopes of 
the ferns may have phylogenetic implication. The basal leptosporangiate fern, 
D. glaucum shares a structure similar to most core-leptosporangiate ferns, but 
differs from that of Pteridaceae sensu Jato indicating that D. glaucum shares 
this feature with the core-leptosporangiate ferns examined to date except 
Pteridaceae sensu Jato. 


Nuclear Evaginations and the Cytoplasmic Organelles 

The present investigation demonstrates that D. glaucum forms nuclear 
evaginations in the mature egg stage of oogenesis. From previous observations, 
it has been shown that Lygodium japonicum also produces nuclear 
evaginations during oogenesis (Cao et al., 2017). However, the fern Osmunda 
japonica (Cao, Dai, and Wang, 2012a), a clade at the base of the 
leptosporangiate ferns, produces no nuclear evaginations during oogenesis. 
Almost all core-leptosporangiate ferns also produce nuclear evaginations 
during oogenesis, including Adiantum flabellulatum, Plagiogyria euphlebia, 
Pteridium aquilinum, Coniogramme emeiensis, Cibotium barometz, Phymato- 
sorus hainanensis, and Anisocampium sheareri (Cao et al. 2010b; Cao, Dai, 
and Wang, 2011; Cao, Dai, and Wang, 2012b; Cao et al., 2017; Wang et al., 
2012a, b; Yang et al., 2013; Zou et al., 2014). Nuclear evaginations are 
considered to be relevant to the information exchange between nucleus and 
cytoplasm, and are important for the latter sporophytic growth (Bell and 
Duckett, 1976). However, the ferns (Osmunda and Equisetum) without nuclear 
evaginations during oogenesis can produce sporophytes as well (Cao, Dai, and 
Wang, 2012a; 2015). Therefore, the biological significance of forming nuclear 
evaginations still needs to be explored. 

The cytoplasmic organelles of D. glaucum are basically similar to those of 
the core leptosporangiate ferns (Cao, Yang, and Wang, 2009; Cao, Wang, and 
Bao, 2010a; Cao et al., 2010b; Cao, Dai, and Wang, 2011; 2012b; Cao et al., 
2017; Wang et al. 2012a, b; Yang et al. 2013; Zou et al. 2013). However, 
vacuolation of the cytoplasm, especially in the upper part of the egg, is 
scarcely seen in other leptosporangiate ferns. 
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The Characteristics of Oogenesis and Their Phylogenetic Significance 
Current hypotheses of fern phylogeny are based on a combination of 
morphological evidence and nucleic acid sequence data (PPG I, 2016; Shen et 
al., 2017; Smith et al. 2006). The morphological evidence of some sporophytic 
reproductive structures (such as the annulus and stomium) have been 
considered to be among the most conserved and consistent characters in 
classification and phylogeny (Shen et al., 2017; Smith, 1995). For the most part, 
anatomical aspects of gametophytic reproductive structures have not been 
included in constructing phylogenetic hypotheses primarily because of the 
technical challenges involved in observing cellular details. Our previous 
observations from the basal groups to the higher groups of leptosporangiate 
ferns have demonstrated that the characteristics of fern oogenesis are similarly 
conserved within orders, but vary between orders (Cao, Yang, and Wang, 2009, 
Cao, Wang, and Bao, 2010a, Cao et al., 2010b, Cao, Dai, and Wang, 2011, 2012b; 
Cao et al. 2017; Wang et al. 2012a, b; Yang et al. 2013; Zou et al. 2014). Thus, the 
evidence to date indicates that fern oogenesis may contribute to understanding 
phylogenetic relationships among clades. The gleichenioid fern D. glaucum, 
which is considered ancestral to the schizaeoid ferns according to the 
phylogenies based on molecular data (PPG I, 2016; Shen et al., 2017; Smith et 
al., 2006), has formed a typical egg envelope and the fertilization pore like the 
core-leptosporangiate lineage. Furthermore, the similar structure of the egg 
envelope and fertilization pore of D. glaucum and most members of the core- 
leptosporangiate ferns except members of the Pteridaceae sensu lato may 
indicate a closer relationship between gleichenioid fern and the core- 
leptosporangiate ferns (cf., Bierhorst, 1971; Holttum, 1973; Lovis, 1977) or 
may show that these features developed independently in the two lineages. 
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ABsTRACT.—Alpine vegetation is one of the most restricted ecosystems in Mexico. Its greatest extent 
occurs on the high peaks of the Trans-Mexican Volcanic Belt (TMVB). In this study we document 
fern diversity of the alpine vegetation of central Mexico, generate a floristic list, provide a key for 
determination, and report the elevational limits of the species present. The latter information will 
serve as a baseline for future comparisons to determine potential upward migration of fern species, 
as is predicted with global warming. In comparison to the paramo vegetation of Central and South 
America, fern diversity in the alpine vegetation of central Mexico is surprisingly low, and only five 
families, 11 genera, and 12 species are documented. This represents 1.5% of the ferns species 
known from the country. Polypodiaceae is the most diverse family, with four genera and four 
species, followed by Pteridaceae (three genera and three species) and Dryopteridaceae (two genera 
and two species). Aspleniaceae and Woodsiaceae each have a single genus, the former with two 
species and the latter with one species. All of the alpine ferns growing in the TMVB are 
leptosporangiate, and there are no arborescent, epiphytic, aquatic, or heliophilous ferns. 
Asplenium castaneum, Cystopteris fragilis, and Polystichum speciosissimum are widespread in 
the alpine vegetation, whereas the remaining species are restricted to one or few sites. With regard 
to upper elevational limit, the highest occurring species were Asplenium castaneum (4,569 m), 
Polystichum speciosissimum (4,490 m), and Cystopteris fragilis (4,377 m). 


Key Worps.—fern diversity, global warming, Mexico, Neotropics, Polypodiaceae 


There are approximately 10,578 species of ferns in the world (PPG I, 2016), 
and nearly a third of them (3,250 spp.) occur in the New World (Moran, 2008). 
Fern communities are richest in wet tropical regions, particularly in cloud 
forests (Kessler, 2010). Mexico is one of the most diverse countries in the world 
for fern species, and Mickel and Smith (2004) report 902 species from the 
country, whereas Villasenor (2016) lists 932 species. Since this latter 
publication, at least one additional species has been described (Sundue, 
2017). Together with Brazil and the Andean region of South America, Mexico 
is considered to represent one of the major centers of fern diversification in the 
Americas (Tryon, 1972). Ferns grow in all of the types of vegetation present in 
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Mexico (Rzedowski, 2006). The greatest richness has been registered in cloud 
forests (approx. 630 spp.), followed by tropical deciduous forests (536 spp.), 
and then tropical rainforests (438 spp.) (Tejero-Diez, Torrez-Diaz and Gual- 
Diaz, 2014). They are also distributed in structurally less complex ecosystems 
such as grasslands; however, they are much less diverse and a mere 21 species 
have been reported (Tejero-Diez, Torrez-Diaz and Gual-Diaz, 2014). Most of the 
diversity is in mountainous areas with a tropical climate, where they occur 
across broad elevational and climatic ranges. 

Alpine vegetation is one of the most restricted ecosystems in Mexico. It 
occurs above timberline at the summits of the highest and coldest mountains. 
This is usually + 3,900 m, but can be as low as 3,600 m in northeastern Mexico 
and on El Ajusco. According to Rzedowski (2006), the Mexican alpine regions 
have mean annual temperatures from 3 to 5°C and can experience frosts 
throughout the year. Average rainfall varies from 600 to 800 mm, and high 
insolation coupled with strong winds result in high evaporation. Alpine 
vegetation occurs in Chiapas (Miranda, 1952) and various locations along the 
Sierra Madre Oriental in the states of Coahuila, Nuevo Leén y Tamaulipas 
(Beaman and Andresen, 1966; McDonald, 1990). It has also been reported for 
Oaxaca (McDonald, 2013) and Chihuahua (McDonald, Martinez, and Nesom, 
2011), but these latter two sites are probably better considered as subalpine. 
Certainly, the greatest extent of alpine vegetation in Mexico is found in the 
Trans-Mexican Volcanic Belt (TMVB), an active volcanic zone that extends 
east-west across central Mexico from the state of Veracruz to the states of 
Jalisco and Nayarit (Gomez-Tuena, Orozco-Esquivel, and Ferrari, 2007). 
Within this region, alpine vegetation is restricted to 11 of the highest 
mountains: Nevado de Colima (NC), Nevado de Toluca (NT), El Ajusco (AJ), 
Iztaccihuatl (IZ), Popocatépetl (PP), Cerro Telapon (TE), Monte Tlaloc (TL), La 
Malinche (LM), Pico de Orizaba, including Cerro Chichimeco (PO), Sierra 
Negra (SN), and Cofre de Perote (CP) (Steinmann et a/., unpublished data). 

Eleven species of ferns have been previously reported from the alpine 
grasslands of the TMVB (Almeida-Lenero et al., 2007; Arreguin-Sdnchez, 
Fernandez, and Quiroz, 2004; Matuda, 1956; McVaugh, 1992; Mickel and 
Smith, 2004; Narave, 1985). However, there is no synthesis of fern diversity 
across the region, and some of the mountaintops remain poorly explored (e.g., 
Sierra Negra and Tlaloc). Therefore, the purpose of this study is to document 
fern diversity of the alpine vegetation of central Mexico, generate a floristic list, 
provide a key to their determination, and report the elevational limits of the 
species present. The latter information will provide a baseline to allow future 
comparisons and determine potential upward migration of fern species, as is 
predicted to occur with global warming. 


MATERIALS AND METHODS 


Bimonthly field trips to the highest mountains of the TMVB were conducted 
over the course of two years (September 2012 to August 2014); sporadic visits 
have occurred since then. We followed Beaman’s (1962) concept of timberline, 
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and only open areas above timberline are included in our study; grasslands 
occurring at lower elevations in openings within otherwise continuous forest 
were not considered. The preparation and collection of herbarium specimens 
followed the techniques presented by Lot and Chiang (1986). The specimens 
were generally gathered from elevations of 3,900 m and higher. Longitude and 
latitude coordinates, as well as altitude, were obtained from a GPS Garmin 
eTrex-30. In addition to the fieldwork, the collections of the following herbaria 
were examined: ENCB, FMCE, IEB, MEXU, MSC, NY, UAMIZ, and XAL 
(Thiers, cont. updated). A dichotomous key for the determination of the TMVB 
alpine fern species was prepared using diagnostic features visible on 
herbarium specimens. Classification and species circumscription follow 
Mickel and Smith (2004), Sigel et al. (2014), and PPG I (2016). Voucher 
specimens collected during the study were deposited in IEB. 


RESULTS 


Five families, 11 genera, and 12 species are here reported from the alpine 
grasslands of the TMVB (Table 1). Nine of these species were found during this 
study. Only Alansmia spathulata (A.R. Sm.) Moguel & M. Kessler, Asplenium 
fibrillosum Pringle & Davenp., and Melpomene peruviana (Desv.) A.R. Sm. & 
R.C. Moran were not recollected, but they are known by previous collections 
from the Cofre de Perote or Iztaccihuatl, at elevations of 3,950-4,150 m. 
Polypodiaceae is the most diverse family, with four genera and four species. 
This is followed by Pteridaceae (three genera and three species) and 
Dryopteridaceae (two genera and two species). Aspleniaceae and Woodsiaceae 
each have a single genus, the former with two species and the latter with one 
species. There are three new records for the alpine vegetation of the TMVB: 
Pleopeltis polylepis (Roem. ex Kunze) T. Moore from the ecotone with Pinus 
hartwegii Lindl. forests on Popocatépet! (4,089 m) and Myriopteris lendigera 
(Cav.) J. Sm. and Pellaea ternifolia (Cav.) Link from the Sierra Negra (4,237 m) 
where they grew together under the shade of large rocks in an area with 
scattered individuals of Pinus hartwegii. 

Of the fern species previously reported from the alpine vegetation of the 
TMVB, Elaphoglossum mathewsii (Fée) T. Moore is probably the taxon that 
was mentioned as simply Elaphoglossum sp. for the alpine vegetation of 
Iztaccihuatl by Almeida-Lenero et al. (2007). Although we did not encounter 
this species during our fieldwork on Iztaccihuatl, it was previously 
documented from there, as well as the Cofre de Perote and Sierra Negra. 
Two additional species included in Mickel and Smith (2004) have upper 
elevational limits that suggest they occur in alpine vegetation: Dryopteris 
pseudofilix-mas (Fée) Rothm. (3,950 m) and Melpomene moniliformis (Lag. ex 
Sw.) A.R. Sm. & R.C. Moran (4,000 m). However, despite our fieldwork and an 
extensive search in herbaria, we have been unable to locate any specimens that 
confirm their presence in the alpine grasslands. Polypodium californicum 
Kaulf. was reported from the alpine vegetation of the Cofre de Perote (Mickel 
and Smith 2004), but this report was corrected to Polypodium calirhiza S.A. 
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TaBLe 1. Diversity, geographical range and sampling dates of ferns of the FVTM alpine belt, * = 
new records. Site abbreviations: AJ, El] Ajusco; CF, Cofre de Perote; IZ, Iztaccihuatl; LM, La 
Malinche; NC, Nevado de Colima; NT, Nevado de Toluca; PO, Pico de Orizaba; PP, Popocatépetl; 
TE, Cerro Telapon; TL, Monte Tlaloc; SN, Sierra Negra. Country abbreviations: Arg, Argentina; Bol, 
Bolivia; Braz, Brazil; Col, Colombia; CR, Costa Rica; Ec, Ecuador; ES, El Salvador; GA, Greater 
Antilles; Guat, Guatemala; Hisp, Hispaniola; Hond, Honduras; JF, Juan Fernandez Islands; Mex, 
Mexico; Nic, Nicaragua; Pan, Panama; USA, United States of America; Ven, Venezuela. 


FVTM sampling location Geographical 
Family Species and elevation (m) distribution 
Aspleniaceae Asplenium castaneum AJ 3,750; CF 3,930-4,040; Central and 
IZ 3,900-4,200; LM southern. Mex; 
4,082; NC 4,000-4,100; Guat, CR, Pan, 
NT 4,175-4,300; PO Col, Ven, Ec, 
3,950-4,377; PP 3,850- Peru, Bol 
4,082; SN 4,271-4,569; 
TL 4,000-4,034 
Asplenium fibrillosum IZ 3,900 Central Mex 
Dryopteridaceae Elaphoglossum CF 3,930; IZ 3,800-4,100; Central and 
mathewsii SN 4,282 m southern Mex; 
Guat, CR, Col, 
Ven, Ec, Peru, 
Bol, Chile 
Polystichum AJ 3,745; CF 3,930-4,267; Central and 
speciosissimum IZ 3,940-3,950; LM southern Mex; 
3,929; NT 4,140; SN Guat, CR, Pan 
4,370; PO 4,086-4,490; 
PP 3,963 | 
Polypodiaceae Alansmia spathulata CF 3,950; IZ 3,990-4,150 Central and 
southern Mex; 
: Guat 
Melpomene peruviana CF 3,930; IZ 4,100; TL Central and 
4,000 southern Mex; 
Guat, Hond, CR, 
Pan, Col, Ven, 
Sur, Ec, Peru, 
Bol, Braz 
*Pleopeltis polylepis PP 4,089 Northern and 
central Mex; 
USA, Guat 
Polypodium calirhiza CF 3,930-4,025 Central Mex; USA 
Pteridaceae *Myriopteris lendigera SN 4,237 USA, Mex, Guat, 
Hond, Nic, CR, 
Pan, Hisp, Col, 
Ven, Ec 
Gaga marginata AJ 3,660 Mex, Guat, Hond, 
ES, GR, Pan, Col, 
Ven, Ec, Peru, 
Bol, Arg 
*Pellaea ternifolia SN 4,237 USA, Mex, Central 
and South 


America, Hawaii 
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TABLE 1. Continued. 


FVTM sampling location Geographical 
Family Species and elevation (m) distribution 
Woodsiaceae Cystopteris fragilis CF 3,930-4,171; IZ 3,800- USA, Mex, Guat, 
4,200; LM 4,082; NC Hond, ES, Nic, 
3,990; NT 4,020; PO CR, Pan, GA, 
3,950-4,377; PP 3,962- Col, Ven, Ec, 
4,082; TE 4,079; TL Peru, Braz, Bol, 
4,000-4,034; SN 4,282 Chile, JF, Arg, 


Europe, Asia, 
Africa, Hawaii 


Whitmore & A.R. Sm. (Sigel et al., 2014). Also, Lehnert (2007) treats the 
Mexican alpine plants of Melpomene, as M. peruviana, instead of M. 
pilossisima (M. Mart. & Gal.) A.R. Sm. & R.C. Moran as has been previously 
reported. 

Asplenium castaneum Schltdl. & Cham. occurred in 10 of the 11 sites 
surveyed, Cystopteris fragilis (L.) Bernh. in nine sites, Polystichum specio- 
sissimum (A. Braun ex Kunze) Copel. in eight sites, Elaphoglossum mathewsii 
and Melpomene peruviana in three, Alansmia spathulata in two, and 
_ Asplenium fibrillosum, Myriopteris lendigera, Gaga marginata (H.B.K.) F.W. 
Li & Windham, Pellaea ternifolia, Pleopeltis polylepis, and Polypodium 
calirhiza were found in a single site (Table 1). The alpine areas with the 
highest fern diversity were the Cofre de Perote and Iztaccihuatl, each with 
seven species, and the Sierra Negra with six species. These were followed by 
El Ajusco, La Malinche, Monte Tlaloc, Nevado de Toluca, and the Pico de 
Orizaba with three species. Popocatépetl and the Nevado de Colima had two 
species, whereas Cerro Telapon was the least diverse with a single species. 
With regard to upper elevational limit, the highest occurring species were 
Asplenium castaneum (4,569 m), Polystichum speciosissimum (4,490 m), and 
Cystopteris fragilis (4,377 m), and all of these limits were reached along the 
mountainous chain of the Pico de Orizaba-Sierra Negra in the eastern portion 
of the TMVB (Table 1). 

All fern species registered by us grow in the cracks of rocky outcrops mainly 
in moist micro-habitats. Only Asplenium castaneum and Cystopteris fragilis 
(Fig. 1) occur in moist microhabitats of protected areas between rocks as well 
as dry sites. They are also the most widely distributed species of the alpine 
grassland and occur both in low and high alpine zones. Finally, we collected 
six species in the low alpine zone: Pleopeltis polylepis was found in the 
ecotone with Pinus hartwegii forest at Popocatépetl, whereas Myriopteris 
lendigera, Elaphoglossum mathewsii, Pellaea ternifolia were encountered at 
the Sierra Negra and Polypodium calirhiza at Cofre de Perote growing under 
the shade of large rocks in an area with scattered individuals of Pinus 
hartwegii. 
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Fic. 1. Habit and close up of the fronds of Asplenium castaneum and Cystopteris fragilis. A. 
Habitat of Asplenium castaneum. B. Close up of Asplenium castaneum. C. Cystopteris fragilis. D. 
Close up of Cystopteris fragilis (Photographs by Victor Steinmann). 


DISCUSSION 


Depending on the origin and features of the sporangium, ferns are divided 
into two groups: eusporangiate and leptosporangiate. The leptosporangiate 
ferns reach their greatest diversity in tropical regions but also are present in 
temperate forests and arid zones. They can be terrestrial, epiphytic, hemi- 
epiphytic, aquatic, or subaquatic. Their size varies from diminutive plants a 
few centimeters tall, such as species of Azolla Lam., to arborescent forms of 
Cyathea Sm. and Dicksonia L’Hér. that exceed 10 meters tall (Mendoza-Ruiz 
and Ceja-Romero, 2014). All of the alpine ferns growing in the TMVB are 
leptosporangiate, low statured, and with small pinnae. The families present 
are among the 10 most diverse fern families in Mexico, in particular the three 
most species-rich families Pteridaceae, Dryopteridaceae, and Polypodiaceae 
(Mickel and Smith, 2004). Although Aspleniaceae and Woodsiaceae are the 
least diverse families, each with a single species, their representatives, 
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Asplenium castaneum and Cystopteris fragilis, respectively, are the most 
widely distributed and abundant ferns in the alpine grasslands surveyed. No 
epiphytic, arborescent, or aquatic ferns occur in the area, and in general, the 
alpine ferns of the TMVB grow in the cracks of isolated, rocky outcrops. 

Only one alpine fern species is endemic to Mexico (Asplenium fibrillosum), 
whereas in total there are about 190 fern species endemic to the country 
(Mickel and Smith, 2004). Thus, the proportion of endemism in the alpine 
vegetation regions is lower than the national average (8% vs. 20%). None of the 
species reported here are strictly alpine, and they all descend to at least 3,000 
m, with some occurring at elevations of less than 1,500 m (Mickel and Smith, 
2004), see Fig. 2. All of the non-endemic species also occur in alpine areas 
outside of the country. In Mexico, Asplenium castaneum has the highest 
distribution; it was found at 4,569 m on the Pico de Orizaba. This species was 
previously known from up to 4,000 m in the region of Iztaccihuatl- 
Popocatépetl and nearly 4,200 m at the summit of the Nevado de Colima 
(Almeida-Lenero et al., 2007; Arreguin-Sanchez, Fernandez, and Quiroz, 2004; 
McVaugh, 1992). Cystopteris fragilis has been found at 4,377 m on the Pico de 
Orizaba. This species is noteworthy for holding the world record for the 
highest occurring fern and has been found growing at 5,485 m in Chang La, 
Ladakh, India (Webster, 1961). It is also considered one of the most 
phenotypically plastic and wide-ranging species in the world, and it occurs 
in numerous mountainous ecosystems with a wide range of altitudes and 
temperatures (Kessler, 2010). These three species are also the most abundant 
and widespread species in the TMVB. 

The alpine vegetation of Central and northern South America is commonly 
called paramo, and in contrast to Mexico, South America has vast expanses of 
alpine vegetation and timberline commonly occurs as low as 3,000 m; it is also 
more humid (Luteyn, 1999). These paramos by far house the greatest diversity 
of alpine plants in the world (Barrington, 2005; Luteyn, 1999), and alpine ferns 
similarly reach their highest diversity there (Tryon and Tryon, 1982), with 
Sklenar et al. (2005) reporting 20 families, 50 genera, and 345 species. In 
comparison, the alpine vegetation of Mexico is poor in ferns, and of its 
approximately 930 species, only approximately 1.5% occur in alpine habitats. 
Seven species are shared between the paramos of South America and the 
alpine regions of the TMVB (Table 1). One of these, Elaphoglossum mathewsii, 
is noteworthy for occurring up to 5,100 m in South America (Luteyn, 1999), 
but in Mexico it occurs to 4,282 m. Melpomene peruviana also reaches a 
considerably higher elevation in South America, 5,200 m as compared to 4,100 
m (Lehnert, 2007). This species is also remarkable because it holds the world’s 
record for the highest elevation vascular epiphyte, recorded at 4,542 m 
(Sylvester et al., 2014). However, in Mexico it is not known to be epiphytic 
(Lehnert, 2007). 

Eleven species are shared with the highlands of Central America (Table 1). 
In contrast, only three species extend northward into the United States. Two of 
the latter are very widespread and also occur in the paramos of South America. 
Thus, the affinities of the alpine fern flora of México are clearly with the alpine 
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Fic. 2. Comparison of elevational ranges of the fern species under study. Ac, Asplenium 
castaneum; Af, Asplenium fibrillosum; Em, Elaphoglossum mathewsii, Ps, Polystichum 
speciosissimum; As, Alansmia spathulata; Mp, Melpomene peruviana; Pp, Pleopeltis polylepis; 
MI Myriopteris lendigera; Gm, Gaga marginata; Pt, Pellaea ternifolia; Cf, Cystopteris fragilis. 
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neotropical vegetation of Central and South America and not Nearctic. 
However, considering that ferns in general are diverse in tropical regions this 
southerly affinity is not surprising. 

There are no studies focusing on the ferns in other alpine areas of Mexico, 
but judging by reports in the available literature, they are even more poorly 
represented than in the TMVB. Garcia-Arévalo and Gonzalez-Elizondo (1991) 
did not encounter any fern species in their survey of alpine regions of Cerro 
Potosi. Farther to the south, Pérez-Ferrera (pers. comm., 2018) reports only 
three alpine fern species in his checklist of the Tacana Volcano in Chiapas: 
(Asplenium castaneum, Polystichum speciosissimum, and an undetermined 
species of Cheilanthes). 

In contrast with other Neotropical alpine zones, for example the Costa Rican 
paramo (Barrington, 2005), the alpine belt of the TMVB harbors no arborescent, 
epiphytic, aquatic, or heliophilous ferns, and in general, the alpine ferns grow 
in moist, often shady cracks of rocky outcrops. It is interesting to note that both 
A. castaneum and C. fragilis occupy from moist to xeric microhabitats in the 
TMVB alpine belt, so their success in the alpine region is probably due to 
characteristics that allow them to survive xeric conditions. 

Additionally, one question to consider is whether the diversity of ferns in 
the alpine belt of the TMVB is increasing due to the rise of species due to 
global warming, as has been registered for other groups of plants (Steinbauer et 
al., 2018). However, there is not a sufficient number of records of specimens of 
herbarium or literature data that allow comparisons. 

The information presented here will permit future comparisons with the 
current diversity and elevational limits of the Mexican alpine ferns. Such 
studies are important because they can be used to detect potential changes in 
the face of global warming (e.g., Pouteau et al., 2016). However, clearly more 
fieldwork is desirable for the alpine areas outside of the TMVB and there are 
many questions to be addressed within the region. For example, compact 
growth forms and morphological characteristics such as buried rhizomes, 
thick fronds with rounded segments, as well as thick indument formed by 
the scales have been correlated with the frigid environments that occur in 
high tropical mountains (Barrington, 2005; Kluge and Kessler, 2007; Tryon 
and Tryon, 1982). Because most of the Mexican alpine ferns have wide 
elevational distributions, it would be interesting to compare morphological 
variation across various ecosystems in order to determine if plants from 
alpine habitats differ from those growing in less extreme conditions and 
how. Furthermore, it would be beneficial to analyze the diversity of Mexican 
ferns across gradients from lower elevations and latitudes to higher 
elevations and latitudes to compare abundance and diversity on a much 
broader scale. 


KEY TO THE ALPINE FERNS OF THE TRANS-MEXICAN VOLCANIC BELT. 


1. Base of the petiole with two vascular bundles; sori with an indusium ............... 2 
1. Base of the petiole with one, three or more vascular bundles; sori without an indusium .. 7 
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2. Petiole straw-colored to red-brown; fronds pinnate-pinnatifid; indusium cup-shaped or 


OVATG ph hi ne oth wh ep Ssipiduean Sid htc: Susmieyopanen< fed a pues BAe eee eee Cystopteris fragilis 
2. Petiole brown, black-purple or blackish; fronds 1-many times pinnate; indusium 
submarginal, elongate 2.0... 6. e ee ce ee ve eee es Che epee 60 3 
3. Sori submarginal, elongate; indusium entire, fibriate or papillate ................... 4 
3. Sori on veins, elongate; indusium entire to suberose ..... 2... ....00).). 20) ee 6 
4¢-¢Pinnae;termate' ew ewepae Ieee ete pas} I OF a Pellaea ternifolia 
4. Pinnae not termate 6. ie. bon neue wheelie pellel ony ciel epepeungs Suey le + opens dane eer 5 
5. Petiole dark castaneous, terete, with orange hairs ................ Myriopteris lendigera 
5. Petiole castaneous, grooved, glabrous’). 5 a.) a... dee som 6m. oo Gaga marginata 
6. *Indusia conspicuously cilate’or laciniate 22-0) goa Geo nem ee Asplenium fibrillosum 
6. Indusia entire to crenulate or somewhat erose, but not ciliate ...... Asplenium castaneum 
7. Fronds simple to pinnatifid .. 02. ..46. 0065. Ge. «tm» vole bt ba tye 8 
7. Fronds once to various times pinnate... i... 5... a ww ss ol «ane te» 10 
8. Fronds pinnatisect to pinnatifid, not dimorphic; rhizome scales concolorous, brown; 
petiole and rachis plaprous sy 0.5. 4.48. 0 ss ee dies Polypodium calirhiza 
8. Fronds simple, dimorphic; rhizome scales bicolorous, to clathrate or red-brown to black; 
petiole and rachis: pubescent’... 0.0. 00.000. 00.0. oe Sle 2? 2 A 9 


9. Fronds dimorphic; rhizome scales red-brown to black; sori acrostichoid ............ 
fF ee eee ON Pe OC Te ee Pe Marge. eae, SE GT Elaphoglossum mathewsii 
9. Fronds not dimorphic; rhizome scales bicolorous to clathrate; sori rounded.......... 
ees es). Saree een ane eit, Fiesty Pleopeltis polylepis 
10. Fertile portion of the frond 2-pinnate to 2-pinnnate-pinnatifid; rhizome scales 
concolorous; petiole and rachis with scales; sori with a false indusium, arranged 


submarginally suite. soma ww hs vet beotetet eee Polystichum speciosissimum 
10. Fertile portion of the frond lobed to pinnate; rhizome scales bicolorous to clathrate and 
iridescent; petiole and rachis with setae; sori otherwise ......... 05 ).90m00 ee 11 


11. Rhizome scales conspicuously clathrate; sori 4-7 per pinna; sporangia with setae..... 
Sn CE Oe PT oe tats. ae oe oft Melpomene peruviana 

11. Rhizomes scales not or scarcely clathrate; sori 1-3 per pinna; sporangia glabrous. ..... ' 
aoeuides Desigialana éoestepage bees flees aed een Alansmia spathulata 


EXSICCATE. 
Alansmia spathulata (A.R. Sm.) Moguel & M. Kessler 


SPECIMENS EXAMINED.—COrRE DE PErote: Veracruz, Mpio. Perote, 3,950 m, Narave 
823 (MEXU). IzTACCiHUATL: Puebla, south side of Iztaccfhuatl, ca. 7 km north of 
Paso de Cortes, 4,150 m, Beaman 2876 (MSC). 


Asplenium castaneum Schltdl. & Cham. 


SPECIMENS EXAMINED.—COFRE DE Perot: Veracruz, Mpio. Perote, ladera suroeste 
del volcan Cofre de Perote (19°29'31.4’”N, 97°08'59.4”O), 4,123 m, Herndndez- 
Cardenas & Arredondo-Amezcua 1526 (IEB). E: Ayusco: Ciudad de México, 
Parque Nacional Cumbres del Ajusco, cordillera oeste, ruta que sube desde el 
albergue alpino (19°12'52”N, 99°15'54”0) 3,750 m, Rincon & Morales 4011 (IEB). 
IzracciHUATL: Estado de México, south side of Iztaccfhuatl, 4,000-4,500 m, 
Beaman 1962 (MSC). La Matincue: Tlaxcala, Mpio. Huamantla, ladera norte del 
volcan Matlalcueye (La Malinche) (19°14’7.05”N, 98°2'7.9”0), 4,082 m, 
Herndndez-Cardenas & Arredondo-Amezcua 1578 (IEB). Monte TLALoc: Estado 
de México, Mpio. Ixtapaluca, ladera S del Monte Tlaloc (19°24'36.5’N, 
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98°42'53.9"0), 4,034 m, Herndndez-Cardenas & Arredondo-Amezcua 1150 
(IEB). Nevapo DE Co.ima: Jalisco, Mpio. Zapotlan el grande (Cd. Guzman), ladera 
noroeste del Nevado de Colima, camino al picacho (19°33'41.3’N, 
103°36'38.1"0), 4,063 m, Herndndez-Cardenas & Arredondo-Amezcua 800 
(IEB). Nevapo DE To.uca: Estado de México, Mpio. Toluca, en la base de la cima 
rocosa del Pico del Aguila, ladera de pendiente pronunciada y de fragmentos 
medianos, cara SO del Nevado de Toluca, 4,300 m, Madrigal & Gonzalez MAGT 
499 (IEB). Pico bE OrizaBa: Veracruz, Mpio. La Perla, ladera sureste del volcan 
Pico de Orizaba (Citlaltépetl) (19°00'08.7”N, 97°16/13.1"0), 4,377 m, Herndn- 
dez-Cardenas et al. 1721 (IEB). PopocarépretL: Estado de México, Mpio. 
Amecameca, ladera norte del Volcan Popocatépetl a medio kil6metro del 
antiguo albergue Tlamacas (19°03’8.6”N, 98°38'7.9”O), 4,082 m, Herndndez- 
Cardenas & Arredondo-Amezcua 836 (IEB). SteRRA NeGra: Puebla, Mpio. 
Atzitzintla, ladera NE del volcan Tliltépetl (Sierra Negra) (18°59/10.4’N, 
97°18'47.6"0), 4,569 m, Herndndez-Cardenas & Arredondo-Amezcua 1014 
(IEB). 


Asplenium fibrillosum Pringle & Davenp. 


SPECIMENS EXAMINED.—IzTAccinuaTL: Estado de México, Mpio. Amecameca, La 
Joya de Alcalican, extremo SW del Iztaccihuatl, 3,900 m, Herndndez 15/78-117 
(IEB). 


Myriopteris lendigera (Cav.) J. Sm. 


SPECIMENS EXAMINED.—SIERRA NgGRA: Puebla, Mpio. Atzitzintla, ladera noreste 
de la Sierra Negra (Tliltépetl), saliéndose a los 2.5 km del camino que lleva al 
GTM, partiendo de la pluma (18°58’55.0"N, 97°18/15.5”O), 4,237 m, Herndn- 
dez-Cardenas & Arredondo-Amezcua 1734 (IEB). 


Cystopteris fragilis (L.) Bernh. 


SPECIMENS EXAMINED.—CERRO TELAPON: Estado de México, Mpio. Ixtapaluca, cima 
del cerro Telapon (19°22'15.6”N, 98°43'12”0), 4,079 m, Hernandez-Cardenas & 
Arredondo-Amezcua 1245 (IEB). Corre DE Perore: Veracruz, Mpio. Perote, ladera 
suroeste del volcan Cofre de Perote (19°29'17.1”N, 97°09’3.7”0), 4,171 m, 
Herndndez-Cardenas & Arredondo-Amezcua 1547 (IEB). IzraccinuatL: Estado de 
México, alpine zone on SW slope of Volcan Iztaccihuatl, 5-6 km N of Paso de 
Cortes, supur of Cerro Amacuilecatl, above and NE of La Joya shelter (ca. 
19°08/15”"N, 98°38'40”W), 4,050 a 4,200 m, Barkley et al. 1038a (MSC). La 
Ma .incue: Tlaxcala, Mpio. Huamantla, ladera norte del volcan Matlalcueye (La 
Malinche) (19°14’7.05”N, 98°2'7.9”O), 4,082 m, Herndndez-Cardenas & Arre- 
dondo-Amezcua 1580 (IEB). Montr TLALoc: Estado de México, Mpio. Ixtapaluca, 
ladera S del Monte Tlaloc (19°24’36.5"N, 98°42'53.9”O), 4,034 m, Herndndez- 
Cardenas & Arredondo-Amezcua 1135 (IEB). Nevapo DE Couima: Jalisco, Mpio. 
Zapotlan el Grande (Cd. Guzman), ladera noroeste del Nevado de Colima, camino 
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al picacho (19°33’47.8’N, 103°36'51.5”0), 3,990 m, Herndndez-Cardenas & 
Arredondo-Amezcua 810 (IEB). Nevapo bE Toiuca: Estado de México, Mpio. 
Tenango del Valle, pendiente sobre barranca, en la ladera E del volcan Nevado de 
Toluca, casi en los limites del bosque, 4,020 m, Eugenia & Héctor MAGT 427 
(IEB). PorpocaréretL: Estado de México, Mpio. Amecameca, ladera norte del 
volcan Popocatépetl a medio kilémetro del antiguo albergue Tlamacas 
(19°03’8.6"N, 98°38'7.9”0), 4,082 m, Herndndez-Cardenas & Arredondo-Amez- 
cua 837 (IEB). Pico pE Orizasa: Puebla, Mpio. Atzitzintla, ladera suroeste del 
volcan Pico de Orizaba (Citlaltépetl) (19°00’27.5"N, 97°17/14.9”0), 4,291 m, 
Herndndez-Cardenas et al. 1472 (IEB). Sierra Necra: Puebla, Mpio. Atzitzintla, 
ladera NE del volcan Tliltépetl (Sierra Negra) (18°58'57.9”N, 97°18’21.9”O), 
4,282 m, Herndndez-Cardenas & Arredondo-Amezcua 1035 (IEB). 


Elaphoglossum mathewsii (Fée) T. Moore 


SPECIMENS EXAMINED.—COFRE DE PEROTE: Veracruz, east side of Cofre de Perote, 
3,930 m, Beaman 2175 (MEXU, MSC). IzraccinuaTL: Estado de México, ladera 
W del Iztaccihuatl, arriba de El] Salto, Valle de Ayoloco, 4,000 a 4,100 m, 
Rzedowski 23469 (MSC). Sierra Necra: Puebla, Mpio. Atzitzintla, ladera 
noreste de la Sierra Negra (Tliltépetl), saliéndose a los 2.5 km del camino que 
lleva al GTM, partiendo de la pluma (18°59/10.6’”N, 97°18/18.9”O), 4,282 m, 
Herndndez-Cardenas & Arredondo-Amezcua 1742 (IEB). 


Gaga marginata (H.B.K.) F.W. Li & Windham 


SPECIMEN EXAMINED.—EL Ayusco: Ciudad de México, Parque Nacional Cumbres 
del Ajusco, cordillera E, ruta que sube de Santo Tomas Ajusco (19°12’46’N, 
99°14'41”0), 3,660 m, Rincon & Cervantes 4072 (UAMIZ). 


Melpomene peruviana (Desv.) A.R. Sm. & R.C. Moran 


SPECIMENS EXAMINED.—COFRE DE PEROTE: Veracruz, east side of Cofre de Perote, 
3,930 m, Beaman 2173 (MSC). IzracciHuaTL: Estado de México, south side of 
Iztaccithuatl, near La Joya, 4,100 m, Beaman 3575 (MSC). Monte TLALoc: Estado 
de México, (19°24’N, 98°43’0), 4,000 m, Wendt & Atkinson (ENCB, NY). 


Pellaea ternifolia (Cav.) Link 

SPECIMENS EXAMINED.—SIERRA NEGRA: Puebla, Mpio. Atzitzintla, ladera noreste 
de la Sierra Negra (Tliltépetl), saliéndose a los 2.5 km del camino que lleva al 
GTM, partiendo de la pluma (18°58'55.0”N, 97°18/15.5”O), 4,237 m, Herndn- 
dez-Cardenas &Arredondo-Amezcua 1733 (IEB). 


Pleopeltis polylepis (Roem. ex Kunze) T. Moore 
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SPECIMEN EXAMINED.—PopocaTEPETL: Estado de México, Mpio. Amecameca, 
ladera norte del Volcan Popocatépetl en el antiguo albergue Tlamacas 
(19°03'30.3”"N, 98°38'7.32”0), 4,089 m, Herndndez-Cardenas & Arredondo- 
Amezcua 841 (IEB). 


Polypodium calirhiza S.A. Whitmore & A.R. Sm. 


SPECIMEN EXAMINED.—COFRE DE PEROTE: Veracruz, Mpio. Perote. Cofre de Perote, 
8.5 km al sur de El Conejo, sobre el camino a las torres (19°29'50’N, 
97°09'06”0) 4,025 m, Ramirez-Amezcua & Steinmann 2011 (IEB). 


Polystichum speciosissimum (A. Braun ex Kunze) Copel. 


SPECIMENS EXAMINED.—COFRE DE PEROTE: Veracruz, north side of Cofre de Perote, 
4,267 m, Stoutamire 2774 (MSC). EL Ayusco: Ciudad de México, Parque Nacional 
Cumbres del Ajusco, cordillera central, ruta que sube al Pico del Aguila 
(19°12'52"N, 99°15/22”0), 3745 m, Rincon & Gomez 4041 (IEB). IzTaccinuatt: 
Estado de México, La Joya de Alcalican extremo SW del Iztaccihuatl, 3,950 m, 
Rzedowski 25398 (MSC). La Matincue: Tlaxcala, Mpio. Huamantla. Ladera norte 
del volcan Matlalcueye (La Malinche) (19°14’25.4”N, 98°01'56.3”O), 3,929 m, 
Hernandez-Cardenas & Arredondo-Amezcua 1114 (IEB). Nevapo bE ToLuca: 
Estado de Mexico, shore of large lake in the crater Nevado de Toluca, 4,140 m, 
Beaman 1884 (MSC). PopocatépretL: Estado de México, Mpio. Amecameca, 
ladera norte del volcan Popocatépetl, Tlamacas (19°03’29.16’N, 98°38'7.08”0), 
3,963 m, Herndndez-Cardenas & Arredondo-Amezcua 2048 (IEB). Pico DE 
OrizaBa: Puebla, Mpio. Atzitzintla, ladera suroeste del volcan Pico de Orizaba 
(Citlaltépetl) (19°00'39.9’N, 97°16'48.1”O), 4,490 m, Herndndez-Cardenas et al. 
1451 (IEB). Srerra Necra: Puebla, Mpio. Atzitzintla, ladera este del volcan 
Thiltépetl (Sierra Negra) (18°59'3.41”N, 97°18/30.2”0), 4,370 m, Hernandez- 
Cardenas & Arredondo-Amezcua 1566 (IEB). 
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Asstract.—The American tree fern Dicksonia sellowiana reaches its northern distribution limit in 
Mexico. This species grows only in the highly fragmented remnants of the montane cloud forest 
(MCF). Given the lack of data on ecological aspects of this endangered tree fern species, considered 
under special protection by Mexican legislation, the objectives of the present research were to 
investigate the structure of D. sellowiana populations and to characterize its habitats. The study 
was carried out in seven populations located at different elevations and latitudes in the MCF of the 
Sierra Madre Oriental in Mexico. In each locality, three sampling plots were established to 
characterize the habitat and to record density, height, stem diameter, number of leaves, and 
individual fertility. The resulting values were analyzed with univariate and multivariate statistics. 
In four of the D. sellowiana populations, individuals of low height predominated; these were 
classified as dynamic populations. However, the density of individuals by unit area was low 
compared with some South American populations, probably because of the high degree of 
fragmentation of the MCF. The variation in population density was mainly explained by the 
vegetation cover. Environmental factors such as elevation and slope orientation, however, 
explained better the variation of the basal area of tree fern populations. These results constitute the 
first reference on the current state of Mexican D. sellowiana populations, which is important for the 
management and conservation of this endangered species, and for monitoring changes in the 
environmental conditions of its natural habitat, the MCF of Mexico. 


Key Worps.—Disturbance, extinction risk, montane cloud forest, population ecology, tree ferns. 


The genus Dicksonia comprises 26 tree fern species (PPG-I, 2016; Noben et 
al., 2017), of which 5 are exclusive to the Neotropics, but only D. sellowiana 
Hook. has a wide distribution (Ulloa-Ulloa et al., 2017a). The latitudinal 
distribution limit of this species in the northern hemisphere is in eastern and 
southern Mexico, in the states of Chiapas, Guerrero, Hidalgo, Oaxaca, Puebla, 
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and Veracruz (Villasenor, 2016). However, it is in Central and South America 
(Argentina, Bolivia, Brazil, Chile, Colombia, Ecuador, Peru, and Venezuela) 
where D. sellowiana populations are the most dense and most widely 
distributed (Ulloa-Ulloa et al., 2017). 

Geographic distribution and population structure are two basic character- 
istics that need to be taken into account to investigate and generate 
conservation and management strategies for plant species at risk (Schmitt 
and Windisch, 2006; Aguraiuja et al., 2008; Mehltreter, 2010). It is therefore 
important to add to quantitative data on the current state of populations and 
habitat characteristics at the local and regional level (Mehltreter and Garcia- 
Franco, 2008; Pérez-Paredes, Sadnchez-Gonzalez, and Tejero-Diez, 2014). 
Mexican populations of D. sellowiana grow exclusively in the montane cloud 
forest (MCF), a type of vegetation with a restricted distribution and a high 
degree of fragmentation (Eleutério and Pérez-Salicrup, 2009; Gual-Diaz and 
Rendon-Correa, 2014), characterized by having high humidity and moderate 
temperature conditions. It is estimated that the MCF currently covers 1% of 
Mexico, but of this area, only 0.4% is primary vegetation, and the remaining 
0.6% is secondary vegetation (Gual-Diaz and Rend6n-Correa, 2014), which is 
the reason why the MCF is at risk of disappearing and urgently requires action 
to be taken for its conservation and management (Arriaga et al., 2000). Given 
this scenario, it is necessary to contribute data on the ecological aspects of 
plant species of the MCF. 

The structure of a plant population can be described through the age, size, or 
stage of the life cycle of the individuals (Aguraiuja, Moora, and Zobel, 2004; 
Weber et al., 2015). In arborescent ferns it is complex or even impossible to 
define the exact age of the individuals, and there is still controversy about the 
degree and type of relationship between fertility and the age of the plants 
(Mehltreter and Garcia-Franco, 2008; Schmitt and Windisch, 2012; Pérez- 
Paredes, Sanchez-Gonzalez, and Tejero-Diez, 2014). Therefore, the stem height 
of individuals (ontogenetic stage of the life cycle) continues to be a criterion 
widely used to describe populations in a single census. This method has 
proven to be very useful in a large number of population studies with different 
arborescent fern species, including D. sellowiana (Mehltreter and Garcia- 
Franco, 2008; Schmitt, Schneider, and Windisch, 2009; Schmitt and Wind- 
isch, 2012; Marques and Kuprek, 2014; Pérez-Paredes, Sanchez-Gonzalez, and 
Tejero-Diez, 2014; Weber et al., 2015; Magna-Silva et al., 2017; Mallmann et 
al., 2018). 

Individuals of D. sellowiana grow slowly and can live up to 200 years 
(Weber et al., 2015); they reach heights of 5-10 m and have fronds up to 5 m 
long (Tryon and Tryon, 1982; Mickel and Smith, 2004). Since it is a perennial 
species with high density and biomass, this raises the environmental 
heterogeneity and favors the coexistence or inhibition of growth of other 
animal and plant species at its natural habitats, giving it a structural and 
functional role in the MCF (Blair et al., 2017; Negrao et al., 2017). Mexican 
populations of D. sellowiana are at risk, mainly due to habitat reduction; its 
trunk is harvested to obtain fiber (“maquique”) used as a substrate for growing 
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Fic. 1. Geographic locations of D. sellowiana populations in the central-south region of the Sierra 
Madre Oriental. 


ornamental plants in the state of Veracruz (Mufiz-Diaz, Mendoza-Ruiz, and 
Pérez-Garcia, 2007; Palacios-Rios, 2007), but in South America it is the tree 
fern species most used for this purpose (Palacios-Rios and Mehltreter, 1999; 
Mallmann et al., 2018). 

Dicksonia sellowiana populations show a wide distribution in the MCF, but 
they develop under precarious conditions in vegetation patches, isolated from 
each other because of fragmentation, which is why they are currently in the 
“special protection” category of the Mexican legislation (Semarnat, 2010) and 
their commercialization is regulated by international agreements (Cites, 2012). 
Given the lack of data on the basic ecological aspects of this species, the main 
objectives of this research were (1) to investigate the population structure of D. 
sellowiana and (2) to characterize its habitat in the Sierra Madre Oriental of 
Mexico. 


MATERIAL AND METHODS 


Study area.—Dicksonia sellowiana populations are distributed in the MCF 
of the Sierra Madre Oriental (SMOR) at elevations ranging from 900 to 2200 m 
(Fig. 1) on slopes of 30°-100°; the dominant climates are Cwa (temperate, dry 
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winter, and hot summer) in Agua Blanca, Medio Monte, Tizapan, Tlanchinol 
(Hidalgo), and Zilacatipan (Veracruz); Cwb (temperate, dry winter, and warm 
summer) in Naupan (Puebla) and Cfb (temperate, no dry season, and warm 
summer) in Xocoyolo (Puebla; Peel, Finlayson, and Macmahon, 2007). The 
average annual temperature at sampling sites ranges between 12 and 23°C, and 
the total annual precipitation between 1000 and 3000 mm, although there is a 
season of low rainfall (the sites have 0-4 months with few rainy days each 
year), which is reflected in foliage loss of some tree elements, but the MCF in 
its primary state maintains its foliage all year. Fog is frequent and atmospheric 
humidity is high, with the consequent decrease in luminosity (Rzedowski, 
2006; Gual-Diaz and Rendon-Correa, 2014). 

Field sampling.—Seven populations of D. sellowiana were analyzed, located 
at different latitudes in the MCF of the SMOR (Fig. 1). According to Hanski and 
Simberloff (1997), a local population is defined as a group of individuals found 
in the same habitat and sharing similar environmental conditions. From this 
definition, the criterion to consider them as populations was based on the 
presence of at least ten individuals in forest patches big enough to establish 
three sampling units under the canopy, each measuring 25 X 25 m (625 m°; 
Alvarez-Ziifiiga et al., 2012; Pérez-Paredes, Sdnchez-Gonzalez, and Tejero- 
Diez, 2014). Schmitt, Schneider, and Windisch (2009) and Alfonso-Moreno et 
al. (2011), suggest a minimum area of 400 m* for the characterization of 
populations of tree fern species, so the size of the sampling plots in the present 
study was adequate. A total of 21 sampling units were established: three in 
each population, at least 20 m apart. 

Density and coverage (basal area) are two important population attributes to 
define the role and ecological status of plant species. Changes occurring in 
both, under different environmental conditions, are of particular interest 
because of the implications for management and conservation of species 
considered at risk (Walker, Mehltreter, and Sharpe, 2010; Negrao et al., 2017). 
In each sampling unit, density (every stem with a frond crown was considered 
to be an individual), stem height, stem diameter at 130 cm (DBH), number of 
mature and fertile individuals, and number of live and senescent leaves were 
recorded (Kingston, Waldren, and Smyth, 2004; Jones et al., 2007). Individuals 
shorter than 1.30 m were measured at the basal part of the frond crown, 
without including the base of the petioles. DBH values were used to calculate 
the basal area of each individual based on the formula zr and the total sum of 
individuals per population (Mateucci and Colma, 1982). The estimation of the 
relationship between stem height, stem diameter, and number of live and 
senescent leaves of the individuals was made with regression analysis; among 
the models generated (linear and non-linear), the one that presented the 
highest determination coefficient value (R*) was chosen (Urban et al., 2010). 

Population structure.—Stem height values were used to classify populations 
as dynamic (Dy), normal/standard (No) or regressive (Re), based on the criteria 
(adapted and simplified) from Aguraiuja, Moora, and Zobel (2004) and Pérez- 
Paredes, Sanchez-Gonzalez, and Tejero-Diez (2014): a dynamic population is 
characterized by the predominance of individuals of small size (in the present 
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study, this included all plants < 1.2 m in height); a normal or standard 
population is defined by the presence of individuals of all size classes, but 
with a predominance of intermediate heights (plants 1.2—-4 m in height); and a 
regressive population is characterized by the dominance of tall individuals 
(>4 m in height) and a limited number of small and/or intermediate-sized 
ones. The morphological values of individuals (except density and basal area) 
in the different populations were contrasted with a one-way variance analysis, 
using Tukey’s test for multiple comparisons of means, following a test of 
normality. | 

As an alternative strategy to define the type of population (dynamic, normal 
or regressive) in a more objective and robust way, the cluster analysis 
technique (CA) was used to identify morphological groups (Galvan-Hernandez 
et al., 2018) based on height, stem diameter, and the number of vegetative 
leaves of individuals (the number of senescent leaves was not included in the 
CA, due to its high correlation with the number of vegetative leaves, r =0.74, 
P<0.001). For this purpose, Ward’s method, considered to be a robust 
algorithm (Cao, Bark, and Williams, 1997; Singh, Hjorleifsson, and 
Stefansson, 2011) was used as a criterion for group clustering. The data were 
then analyzed using an inverse ANOVA to find the number of morphological 
groups that explained the total and intrinsic variation of the sample (the 
attributes of the species itself as a basis for class segregation). The clusters 
obtained were contrasted by k means to describe the significant differences 
among them (Stat Soft Inc., 2010). 

Habitat characterization.—In each of the 21 sampling units, the habitat was 
characterized and the degree of disturbance determined by estimating the 
variables proposed by Pérez-Paredes, Sanchez-Gonzalez, and Tejero-Diez 
(2014), modified for the conditions observed during fieldwork. The following 
variables were estimated: distance between sampling sites and human 
settlements (the values were obtained directly from Google Earth Pro, DIST); 
elevation (m, ELV) obtained with a Global Positioning System Garmin eTrex 
10; terrain slope (degrees, SL), calculated with a Suunto clinometer; 
orientation (degrees, OR), recorded with a Suunto Tandem compass; 
percentage of vegetative ground cover (GVC) and percentage of rock cover 
(RC), using the Braun-Blanquet scale (1979); canopy cover (CC), estimated with 
a convex spherical model A densiometer (Lemmon, 1956); number of roads 
(trails [T], unpaved roads [UR], highways [H]) within each plot; felling or 
clearing (EXT) was estimated using a subjective scale, with values ranging 
from 1-10, divided into three intervals: 1-3 representing the lowest values, 4— 
6 the intermediate levels, and 7-10 the high levels; the presence of landslides 
(LS) was classified according to the area they occupied in m* (Appendix). 

The degree of relationship between population density (or basal area values), 
and variable values obtained by habitat characterization was estimated 
through a non-linear analysis technique (regression trees; RT). RTs are very 
effective and useful for expressing non-linear relationships, since they do not 
require a priori assumptions about any type of behavior between the variables, 
and the interpretation of the results is simple, which is not the case with linear 
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statistical models (Sanchez-Gonzalez and L6pez-Mata, 2005; Kallimanis et al., 
2007; Krzywinski and Altman, 2017). 

Recursive partitions of the population density values (or basal area values) 
from the 21 sampling sites were made using the RT algorithm, resulting in 
increasingly homogeneous. subgroups. At each division, the values of the 
categorical or continuous environmental values used to characterize the 
habitat were selected as threshold values that resulted in the greatest reduction 
in the error sum of squares of population density (or basal area) in the 
sampling sites. The resulting RT structure has the following characteristics: it 
begins at a node (root), and continues through a series of binary divisions 
(stems), determined on the basis of values of the independent variables (used 
to characterize habitat) and ends at a series of terminal nodes (leaves). The 
number that appears at the start of the observations that follow the same route 
from the root to a particular leaf represents the arithmetic mean of the 
population density (or basal area) of a particular subset of D. sellowiana 
populations (Kallimanis et al., 2007). All the analyses were calculated using 
the Statistica version 10 math program (Stat Soft Inc., 2010). 


RESULTS 


Characterization of population structure.—A total of 470 individuals were 
counted (3.6 individuals/100 m7), and population density in sample sites 
varied from 31-92 plants. The highest density of individuals was recorded at 
Tizapan (n=92), followed by Medio Monte and Agua Blanca with 86 and 70, 
respectively (Fig. 2). 

The classification of population types based on individual size categories 
was similar to the classification obtained using cluster analysis criteria. The 
cluster analysis resulted in the identification of three significantly different 
categories or morphological groups in the populations (Table 1, Fig. 3). The D. 
sellowiana populations of Medio Monte, Tizapan, Xocoyolo, and Zilacatipan 
(Group 1) showed a high number of individuals (between 54 and 90%) in an 
early development class of less than 1.2 m in height, so they were classified as 
dynamic (Dy). The populations in Agua Blanca, Naupan, and Tlanchinol 
(Group 2) showed a predominance of intermediate class individuals (between 
1.2 and 4 m height), with 78, 51 and 37%, respectively, so they were classified 
as normal (No). 

The population in the municipality of Tlanchinol showed the highest mean 
DBH (31.48+16.81), height (2.48+1.77), number of green fronds (13.84+8.07) 
and number of senescent fronds (13.45+8.74). The lowest values of DBH 
(13.72+7.88), number of vegetative leaves (5.70+3.23) and number of 
senescent leaves (4.18+3.70) were observed in the Xocoyolo population, and 
the lowest mean height at Tizapan (0.77+1.0; Fig. 4). 

There were significant statistical differences among populations (P< 0.01) in 
DBH, height, and number of green and senescent leaves (Fig. 4, Table 2). A 
highly significant positive correlation (P < 0.001) was found between height 
and diameter (R=0.72, y=(0.17+0.03x)*); height and number of vegetative 
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Fic. 2. Dicksonia sellowiana population types defined based on size structure (heights of 
individuals). 


leaves R=0.65, y=(3.1446.25x)°; height and number of senescent leaves 
(hyn. y=(2.6146.16x)*); diameter and number of green leaves (R=0.75, 
y=(1.5+0.06x)*); and between diameter and number of senescent leaves 
(R=0.71, y{1.05+0.07x)°). | 

Fertility—Among the 470 individuals analyzed, 305 were fertile, which 
amounts to 64% of the total—the number of fertile individuals in each 
population varied from 24-59. The percentages of fertile individuals were: 
Naupan (88.1), Medio Monte (81.4), Tlanchinol (77.6), Zilacatipan (77.4), Agua 
Blanca (74.3), Tizapan (46.7), and Xocoyolo (30.8). At Medio Monte, Naupan, 
Tizapan, Xocoyolo, and Zilacatipan some individuals with heights of 0.13— 
0.50 m (a mean of 18 cm in diameter) with developed sporangia were observed. 
Sterile individuals of more than 2 m in height were found in all populations 
(Fig. 2) except Tizapan, where all individuals > 2 m were fertile. 


Tas_eE 1. Morphological groups defined using cluster analysis. The specific classification criteria 
for Dicksonia sellowiana populations are indicated. Letters indicate statistical significant 
differences (p < 0.001) among groups obtained by k-means analysis (F(2,463) = 589.99; F’(2,468) = 
187.03; Fo 468) = 462.04). DBH: Diameter at Breast Height. N: Number of individuals. 


Morphological group DBH Height Vegetative leaves N 
Small individuals 13.79+8.40° 0.70+0.82* 5.24+2.94° 259 
Medium individuals 28.46+4.02” 2.03+1.04” 14.17+4.42”" 187 


Large individuals 48.99+11.85° 3.94+1.51° 21.20+5.42°" 25 
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Fic. 3. Dendrogram of cluster analysis of D. sellowiana populations showing morphological 
groups identified on the basis of DBH (Diameter at Breast Height), height, and number of vegetative 
leaves. 


Relationship between population density and habitat characteristics.—All 
populations of D. sellowiana found during the field trips were growing within 
the forest patches, but not at the periphery, so it was not necessary to consider 
the edge effect in the sampling design. The results of the regression tree 
indicate that ground vegetation cover (GVC), clearing or extraction (EXT), and 
distance between sampling sites and human settlements (DIS) are the variables 
most related to variations in D. sellowiana population density. The lowest 
number of individuals was found in the plots where GVC was less than 62.6% 
and the clearing rate was higher than 1.5; and the highest mean density of 
individuals occurred in sites with values of GVC above 62.6% and distances to 
human settlements less than 13.5 min (in one sampling plot) or less than 36 
min (in four sampling plots) (Fig. 5). 

Relationship between basal area (biomass) and habitat characteristics.—The 
regression tree indicates that elevation and slope orientation were the two 
variables that best explained variation in basal area of the populations. In the 
sense that the basal area tends to increase with elevation, except for three 
plots, corresponding to the population in the Tlanchinol municipality (1458.9 
m a.s.l.), in which the mean BA value was higher than 1000 m*/ha at elevations 
below 1595 m. The basal area was larger in localities with a northwest 
exposure (Fig. 6). 
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Fic. 4. Average values of attributes of D. sellowiana populations in Agua Blanca (AB), Naupan 
(Na), San Bartolo (SB), Tizapan (Ti), Tlanchinol (Tl), Xocoyolo (Xo) and Zilacatipan (Zi). Different 
letters indicate significant differences. 


DISCUSSION 


Density in populations of D. sellowiana is variable and differs by region: in 
Brazil, Gasper et al. (2011), Weber et al. (2015), Mallmann et al. (2018), Biondi 
et al. (2009), and Marques and Krupek (2014) counted 1.5, 1.5, 13.2, 13.3, and 
13.8 individuals (average values per plot, standardized to 100 m°), respective- 
ly, and in Colombia, Alfonso-Moreno et al. (2011) recorded 14.2. The average 
number of D. sellowiana individuals per plot in SMOR was low (3.6 
individuals /100 m7), compared to those recorded in some South American 
populations, and relatively similar to the values obtained for other species in 
Mexico such as Alsophila firma (Baker) D.S. Conant (8.1 individuals /100 m7) 


Taste 2, ANOVA test of four morphological characteristics of D. sellowiana populations in the 
Sierra Madre Oriental, Mexico.*P<0.001 


Variable Df Sum Sq Mean Sq F value 
Height 6 1357 22.618 16i27% 
DHB 6 18132 3022.1 28. 2* 
Vegetative fronds 6 3799 633.2 19.6* 

Senescent fronds 6 4641 773.4 17.44* 
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Fic. 5. Regression tree of population density in D. sellowiana. Circles and squares represent non- 
terminal and terminal nodes, respectively. Number in each circle or square indicates the mean 
number of individuals in the populations varying according to the branches of the tree towards a 
specific terminal node. Number under each circle or square indicates the sum of squares 
(associated with the arithmetic mean of all samples remaining through the nodes). Numbers in the 
links between nodes represent the decision criteria by which a particular environmental variable 
provides the basis for division. Solid black lines represent the route determining the highest 
population density, and gray lines the route determining the lowest population density. GVC. 
Ground vegetation cover; EXT. Clearing or extraction; DIS. Distance to the closest settlement. 


and Cyathea fulva (M. Martens & Galeotti) Fée: 4.4/100 m* (Pérez-Paredes, 
Sanchez-Gonzalez, and Tejero-Diez, 2014). 

Variation in population density has been associated mainly with the degree 
of habitat disturbance (Marques and Krupek, 2014; Pérez-Paredes, Sanchez- 
Gonzalez, and Tejero-Diez, 2014; Magna-Silva et al., 2017). Consequently, the 
low individual density of D. sellowiana by unit area could be related to the 
increasingly fragmented and reduced distribution of MCF in Mexico, mainly 
caused by deforestation, farming, mining, and urban spread (Sanchez-Ramos 
and Dirzo, 2014). In contrast, other populations of large-sized ferns 
(Lophosoria quadripinnata (J.F. Gmel.) C. Chr.) and tree ferns (Cyathea 
divergens Kunze) are favored and show high density in disturbed environ- 
ments (Eleutério and Pérez-Salicrup, 2009). 

It has been suggested that the variation in population density of D. 
sellowiana is due to the intensity of light reaching the understory (Weber et 
al., 2015), to humidity and edaphic conditions (Alfonso-Moreno et al., 2011; 
Mallmann et al., 2018), to low water availability and low anthropic activity 
(Marques and Krupek, 2014), and to low drought periods (Gasper et al., 2011). 
One of the environmental factors most correlated to high individual density of 


36 AMERICAN FERN JOURNAL: VOLUME 109, NUMBER 1 (2019) 


Basal Area 


67729 
ALT>1403 


60762 
ALT>1595 


ALT>2054 


Fic. 6. Regression tree of basal area in D. sellowiana. As explained for Fic. 5, but here the 
dependent variable is mean population basal area. Vectors with black lines represent the route 
determining the largest basal area, and vectors with gray lines the route determining the lowest 
basal area. ELV. Elevation; NE. Northeast; NW. Northwest; ORI. Orientation. 


D. sellowiana in SMOR was ground vegetation cover (GVC). It is probable that 
this variable is an indicator of local environmental conditions suitable for 
successful development of the species (under some degree of disturbance in 
the forest). Another factor related to population density was the distance 
between sampling localities and human settlements. Localities closest to 
settlements showed a high population density, which could seem contradic- 
tory since they would be expected to be more susceptible to human-caused 
disturbances. However, in Mexico most forests are property of private owners 
or “ejidatarios”, which is why the localities closest to towns are better 
supervised and those farther away are not. A good example is the Tlanchinol 
locality, where the owners of the forest routinely visit their lands, so it is 
difficult for other people to extract plants and cut down trees. 
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Another important factor to consider is the environmental history of the 
localities where D. sellowiana grows. In the Tlanchinol population, for 
example, where the most robust tree ferns are found, natural regeneration of 
the forest was common in the past, but in recent years, people clear or “weed” 
the trails they customarily use. In contrast, in the Naupan locality, forest fires 
are common and unknown people furtively extract entire small individuals 
and half the trunk of more robust individuals. In Xocoyolo the forest was 
fragmented and surrounded by pasture. In this locality craftsmen extract the 
“maquique” or “xaxim” (root fiber) of two species of arborescent ferns 
(Alsophila firma and Cyathea fulva), but not of D. sellowiana, since they claim 
that the trunk of this species does not have the necessary characteristics they 
need for their products. The only use of D. sellowiana observed in the present 
study was for sale as ornamental plant at traditional markets in the Sierra Norte 
in the state of Puebla. 

An analysis of the structure of tree fern populations can be useful for 
recognizing patterns at a regional level (Aguraiuja, Moora, and Zobel, 2004). In 
the case of D. sellowiana in the SMOR, the forest was in different stages of 
regeneration, as mentioned in the previous paragraph, so the structure was 
dynamic in four populations, and normal in three, which is consistent with the 
results observed in other populations of the same species in different 
provinces in Brazil and Colombia, where individuals of lower height were 
the most abundant in places with different type of anthropic disturbance 
(Biondi et al., 2009; Alfonso-Moreno et al., 2011; Gasper et al., 2011; Weber et 
al., 2015). Predominance of seedlings, young, mature or old individuals is a 
product of the interaction between multiple environmental factors and/or 
random events at the local level (selective extraction of plants, grazing, 
opening of roads, natural or induced fires, landslides, dry years, among others) 
in a given time and space (Marques and Krupek, 2014; Pérez-Paredes, Sanchez- 
Gonzalez, and Tejero-Diez, 2014), which is why it is common to find tree fern 
populations with different size or age structures. For example, Marques and 
Krupek (2014) recorded regressive populations (in which the number of 
mature individuals was larger than the number of young ones) of D. sellowiana 
in a fragment of forest in good condition (high humidity, low anthropogenic 
activity) in southern Brazil. 

The predominance of individuals of less than 1.2 m in four of the seven 
populations analyzed in the present study could be considered as a criterion to 
define stability or instability of D. sellowiana populations (Aguraiuja, Moora, 
and Zobel, 2004; Odland, 2007; Aguraiuja et a/., 2008), or as an indicator that 
environmental conditions are suitable for spore germination, gametophyte 
formation, sexual reproduction, and seedling development (Large and 
Braggins, 2004; Sharpe and Mehltreter, 2010; Hernan-Cassini, 2013). 

Under controlled conditions, the spore germination rate of D. sellowiana is 
high (88%) and the spore viability is sufficient to generate a spore bank in the 
forest (Fiilippini, Duz, and Randi, 1999; Schmitz et al., 2006). Consequently, it 
is important to identify the environmental factors most correlated with 
survival in other perhaps more critical stages of its life cycle, such as 
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recruitment and subsequent development of seedlings (Sharpe and Mehltreter, 
2010), in order to make suitable proposals for management and conservation of 
the species (Bernabe, Williams-Linera, and Palacios-Rios, 1999). The predom- 
inance of intermediate-sized individuals in the populations classified as 
normal or stable (Agua Blanca, Naupan, and Tlanchinol) could be related to 
moderate disturbance events, either natural or caused by human activities 
(Pérez-Paredes, SAnchez-Gonzdalez, and Tejero-Diez, 2014) such as “selective 
extraction” of small individuals or parts of robust individuals (as in the 
Naupan and Tlanchinol localities), and/or low availability of micro-environ- 
ments for the establishment of seedlings because of factors such as landslides 
or trail-opening (Bernabe, Williams-Linera, and Palacios-Rios, 1999; Sharpe 
and Mehltreter, 2010). 

Fertility in fern populations tends to vary spatially and temporally (Sharpe 
and Mehltreter, 2010)—in the present study, differences were found in fertility 
percentages between populations, especially in the classes under 0.58 m, so it 
is possible either that physiological development at this height is not enough 
for the plants to produce spores or that local environmental conditions are not 
favorable for fertility (Bernabe, Williams-Linera, and Palacios-Rios, 1999; 
Sharpe and Mehltreter, 2010). According to Magna-Silva et al. (2017), it is 
common in tree ferns to find that small-sized, fertile individuals actually 
originate from mature plants by vegetative reproduction. Field observations 
confirm this claim—some individuals of D. sellowiana were growing as aerial 
stems from recumbent and/or partially buried mature individuals. 

Neumann, Schneider, and Schmitt (2014) suggest that Cyathea corcovaden- 
sis (Raddi) C. Christensen has the potential to increase its population density 
when the proportion of fertile individuals is over 60%. Assuming a similar 
reproductive behavior in other tree fern species, the mean fertility percentage 
in D. sellowiana at a regional level (SMOR) was high (67%). However, at the 
local level, the populations of Tizapan and Xocoyolo showed lower fertility 
percentages, 46.7 and 30.8%, respectively, so they could be considered as 
declining according to Neumann, Schneider, and Schmitt’s proposal (2014). In 
two D. sellowiana populations from southern Brazil, “low” fertility values 
were also recorded; 52%, and 38% (Schmitt, Schneider, and Windisch, 2009; 
Marques and Krupek, 2014, respectively), and, in a population in Colombia, 
only the individuals over 4 m in height (44.4%) were fertile (Ramirez-Valencia, 
Sanin, and Alvarez-Mejia, 2009). However, it would be necessary to carry out 
specific studies of the relationship between fertility and population density to 
determine whether the proposal of Neumann, Schneider, and Schmitt (2014) is 
accurate. 

Intermediate and high basal area values (indirect measure of biomass) of the 

D. sellowiana populations occurred at elevations above 1,450 m, which is 
consistent with their climate and habitat preferences, given that their 
distribution is restricted to the MCF (Mickel and Smith, 2004), unlike other 
tree fern species such as Alsophila firma which tends to be successful at 
elevations lower than 1,470 m in the ecotone between the MCF and the sub- 
evergreen tropical forest in the states of Hidalgo and Puebla, in the SMOR, 
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Mexico (Eleuterio and Pérez-Salicrup, 2009; Pérez-Paredes, Sanchez-Gonzalez, 
and Tejero-Diez, 2014). Elevation and latitude are environmental values that 
do not have a direct effect on plants—in contrast, temperature and humidity, 
which change along latitudinal and elevational gradients, do have a direct 
effect on the distribution of species (Bystriakova, Schneider, and Coomes, 
2011; De Frenne et a/., 2013). The populations of D. sellowiana analyzed here 
are found at the northern limit of the neotropical distribution of the species— 
there are no records of tree ferns in the MCF of the SMOR in more northern 
latitudes, such as in the El Cielo Biosphere Reserve in Tamaulipas. It is 
probable that climate factors that vary negatively with latitude and elevation, 
such as temperature and precipitation (Jardel-Pelaez et al., 2014), interfere 
with the abundance and distribution of tree fern species (Bystriakova, 
Schneider, and Coomes, 2011; De Frenne et al., 2013; Sosa et al., 2016). 

The montane cloud forest is considered a highly vulnerable ecosystem on a 
worldwide level, due to its high degree of fragmentation and sparse 
occurrence. In Mexico, the main problem faced by this ecosystem is related 
to changes in land use to crop and livestock farming (Gual-Dfaz and Rendon- 
Correa, 2014). In this scenario, it is increasingly important to study the 
distribution and structure of plant populations that are considered represen- 
tative of the “typical” environmental conditions for this vegetation type. In the 
particular case of D. sellowiana, it is important to analyze the current state of 
its populations, because it is highly linked to its habitat conditions, which 
means that it could be used in restoration programs as well as for monitoring 
changes in forests over time. 
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Asstract.—The present study describes the anatomy of Anemia villosa, Anemia organensis, and a 
hybrid of the two species (Anemia villosa x Anemia organensis). The aim of the present study was 
to answer the following question: can the Anemia villosa x Anemia organensis hybrid be 
distinguished from the parent plants based on its anatomical traits? The results showed that the 
two species and hybrid can be easily differentiated by their anatomical characteristics and that the 
hybrid anatomy shares features with both parents. In addition, the hybrid differed from the parents 
in two variables: leaf thickness and number of trichome cells. 


Key Worps.—Anemia organensis, Anemia villosa, discriminant analysis, hybridization, stele. 


Hybridization is common in plants, especially ferns (Barrington, 1989). 
Anemia is a fern genus of the family Anemiaceae with 115 species in the 
Neotropical region and 17 in Africa, Madagascar, Réunion, and southern India 
(Mickel, 2016). Many Anemia species are apparently of hybrid origin, and 35 
hybrid combinations have been identified in the Neotropical region, with more 
likely to be recorded (Mickel, 1962; 1982; 2016; Pena, Schwartsburd, and 
Alves-Aratijo, 2017; Rabelo and Schwartsburd, 2016). Brazil is the primary 
center of diversity of Anemia with 66 species, of which 39 are endemic 
(Anemia in Flora do Brasil 2020) and 19 hybrids recorded (Mickel, 2016; Pena, 
Schwartsburd, and Alves-Aratijo, 2017; Rabelo and Schwartsburd, 2016; 
Santos and Araujo, 2007). 

Anemia organensis Rosenst. is an endemic species of the Brazilian Atlantic 
Rain Forest (Anemia in Flora do Brasil, 2020) and only one hybrid with 
Anemia villosa Humb. & Bonpl. ex Willd. has been recorded for the 
municipality of Santa Maria Magdalena, Rio de Janeiro State, Brazil (Mickel, 
1962; 2016). In addition to the hybrid with A. organensis, hybrids of A. villosa 
have been recorded with A. raddiana and A. imbricata, all in Brazil (Mickel, 
2016). 

Morphological data are frequently used to verify the occurrence of hybrids in 
ferns (e. g., Kott and Britton, 1982; Barrington, 1986; Paris and Windham, 1988; 
Haines, 2003), and other groups (Fant and Preston, 2004). Anatomical features 
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can also be useful in hybrid identification (e.g., Webb and Carlquist, 1964; 
Brandham and Cutler, 1978; Baas, 1978). According to Baas (1978), anatomy 
deserves greater attention and could be used to verify the occurrence of 
hybrids. Hybridization research on ferns typically considers only morphology 
or cytology (e.g., Kott and Britton, 1982; Barrington, 1986; Paris and Windham, 
1988) and we found no studies in the literature on the anatomy of fern hybrids. 

Some studies involve hybrids that are artificially produced in laboratories 
and compare their anatomy with naturally occurring specimens. This kind of 
research requires a range of data to accurately identify hybrids (Judd et al., 
1999). Multivariate exploratory techniques are useful in distinguishing 
hybrids and principal component analysis (PCA) and discriminant factor 
analysis (DFA) have been used by a number of researchers (Thébaud and 
Abbot, 1995; Gurgerli, 1997; Levin, 1999; Albaladejo, Aparicio, and Silvestre, 
2004). 

The present study describes the anatomy of Anemia villosa, Anemia 
organenis, and a hybrid of the two species (Anemia villosa x Anemia 
organensis). The aim is to answer the following question: can the hybrid 
(Anemia villosa x Anemia organensis) be distinguished from the parent plants 
based on anatomical variables? And if so, what anatomical characters are 
involved? 


MATERIAL AND METHODS 


Leaf samples of Anemia villosa, Anemia organensis and the Anemia villosa 
x Anemia organensis hybrid were collected along the “Pedra Dubois” trail, in 
the municipality of Santa Maria Madalena (21°45’-22°00’S and 41°41’- 
42°00'W), Rio de Janeiro state, Brazil. 

Five individual parent plants and hybrids were randomly selected (N=15). 
Five fully developed leaves were collected from each and immediately fixed in 
70% formalin-acetic acid-alcohol in the field (FAA) (Berlyn and Miksche, 
1976). Next, the material (0.5 cm* of each leaf) was gradually dehydrated in 
ethanol and embedded in glycol methacrylate resin (Feder and O’Brien, 1968). 
Sections (5 1m) were cut using a rotary microtome and subsequently stained 
with 0.1% Toluidine blue O (O’Brien and McCully, 1965). 

Stomatal density was estimated in dissociated material using Jeffrey and 
Franklyn’s reagent (Jensen, 1962). 

Anatomical descriptions and measurements were performed using an 
Olympus BX50 light microscope and image analysis software (Image-Pro Plus, 
version 3). The images were obtained with a Cool SNAP-Pro digital camera 
attached to the microscope. 

For scanning electron microscopy, leaf fragments were fixed in 2.5% 
glutaraldehyde, 4% formaldehyde in 0.05 M cacodylate buffer. After 
dehydration in acetone, the samples were placed in a critical point Bal-Tec 
Critical Point Dryer CPD 030. The fragments were then affixed with carbon 
tape on their own supports and covered with a thin layer of 20 nm gold (Bal- 
Tec Sputer Coater SCD 050). The samples were then observed and electro- 
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micrographed under the scanning electron microscope ZEISS DSEM 962 and 
EVO 40 at a voltage of 25KV. 

The spores were collected in the field, packed in envelopes and mounted on 
slides with glycerin, according to Large and Braggins (1990). 

The following leaf anatomical parameters were measured for each 
population: thickness of mesophyll, palisade and spongy parenchyma, adaxial 
and abaxial epidermis, and leaf (vein and non-vein regions). The sample size 
established was 25 fields for each measurement (N=25), from which the means 
and standard deviations were calsulated. The results of the measurements 
were used in a simple correlation between anatomical traits, and the non- 
anatomical factors were assessed. The significance of differences between the 
populations studied was established by multivariate analysis of variance 
(ANOVA). Principal components analysis (PCA) was used to group the 
populations and identify the components with the highest variance. 
Discriminant Function Analysis (DFA) was applied to verify the hypothesis 
that the study sample consisted of discrete groups that are morphologically 
different from each other (Schlichting, 1986; Zar, 1996). Overall differences 
between the groups are expressed as the Mahalanobis distance. The statistical 
analyses were performed using Statistica v. 6.0 software (Statsoft). 


RESULTS 


The laminae in the parent plants and hybrid are 1-pinnate-pinnatifid in A. 
villosa, 1-pinnate in A. organensis, and a shallowly pinnatifid laminae in the 
hybrid. The spores of Anemia villosa and Anemia organensis are trilete, 
averaging 84.18 tm and 92.02 pm long, respectively. The hybrid spores are 
deformed, with an average length of 81.73. 

The petiole (Fig. 1) of the parent species and hybrid exhibit a uniseriate 
epidermis with lignified walls. The stele is “V” shaped in A. villosa (Fig. 1 A) 
and “C” shaped in A. organensis and the hybrid (Fig. 1 B and C). 

Both species and the hybrid displayed hypostomatic leaves (Fig. 2 B, D, and 
E). A single-layer epidermis was observed, with sinuous anticlinal cell walls 
(Fig. 2 B, D, and E). Stegmata are present on the epidermis cells above the veins 
(Fig. 3 G, H, and J). In the sections, the external periclinal walls of the 
epidermal cells are convex. Multicellular trichomes are present in both species 
and the hybrid. 

The mesophyll in both species and the hybrid contains 3-4 layers of 
branched parenchyma (Fig. 3). The first layer below the adaxial epidermis is 
more compact. In A. organensis (Fig. 3 L) and the hybrid (Fig. 3 K), the first and 
second layers below the adaxial epidermis display elongated cells (Fig. 3 K 
and L). 

The leaf margin differs between the studied species and is straighter in A. 
villosa (Fig. 3 J), but angled toward the. abaxial surface in the hybrid and A. 
organensis. In the hybrid, the leaf margin has a vein consisting only of xylem 
cells (Fig. 3 K). 
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Fic. 1. Petiole, transverse section. (A, D) Anemia villosa (B, E), hybrid, and (C, F) Anemia 
_organensis Bars= 350um (A, B, and C) and 100 pm (D, E and F). E= endodermis, F= phloem and X= 
xylem. 


Table 1 shows the ANOVA results of the groups for the quantitative 
variables. The parent groups differed for 10 of the 11 variables studied. The 
hybrid differed from the parents in two variables: leaf thickness and number of 
trichome cells. The hybrid is similar to A. organensis in five variables (adaxial 
surface epidermal thickness, stomatal length, leaf thickness in the non-vein 
region, mesophyll thickness, and trichome density on the adaxial surface) and 
to A. villosa in four variables (abaxial surface epidermal thickness, stomatal 
density, trichome density on the abaxial surface, and spore length) 

Factor loadings and eigenvalues for the first two components (PCs1 and 2) 
extracted in PCA are shown in Table 2. These accounted for 62% of the total 
variance. Anatomical traits such as leaf thickness, mesophyll thickness, and 
trichome density on the adaxial surface showed the highest (either positive or 
negative) correlations with PC1, and trichome density on the abaxial surface 
and spore length exhibited the highest correlations with PC2. The scatterplot 
(Fig. 4) for these two components indicated three distinct groups. 

Factor structure for the 11 variables and eigenvalues for the first two factors 
(DF1 and 2) extracted in DFA are shown in table 2. DF1 was positively 
correlated with leaf thickness and negatively with stomatal density, whereas 
DF2 was positively correlated with the stomatal density. The DFA classifica- 
tion was 100% for individuals from groups I, II, and III. 
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Fic. 2.. Epidermis, scanning electron micrographs. Anemia villosa on the adaxial (A) and abaxial 
surface (B). Hybrid, on the adaxial (C) and abaxial surface (D). Anemia organensis on the adaxial 
(E) and abaxial surface (F). Bars= 200 um (A, C, E); 20 um (B, D, and F). T= trichomes, — sinuous 
anticlinal wall. 


DISCUSSION 


The results showed differences between the hybrid and parent plants in 
terms of quantitative and qualitative anatomical characters. Mickel (1962) 
described the Anemia organensis x Anemia villosa hybrid considering only the 
shallowly pinnatifid laminae and aborted spores. These morphological 
characteristics are intermediate between the pinnate-pinnatifid laminae in A. 
villosa and 1-pinnate laminae in S. organensis (Mickel, 2016). Moreover, the 


hybrid spores were deformed, whereas A. villosa and A. organensis displayed 
trilete spores. 
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Fic. 3. Main vein region, secondary veins, stegmata and leaf margin, transverse sections (A, D, G, 
and J) Anemia villosa, (B, E, H, and K), hybrid (C, F, I, and L), and Anemia organenis Link. Bars= 
100 um (A, B, C, D, E, and F), 25 pm (G, H, and I) and 50 pm {J, K, and L). — stegmata, VB = vascular 
bundle, T= trichome, CPW= convex periclinal wall. 


Baas (1978) studied hybridization in I/ex species and identified four possible 
hybrid anatomies: 1) intermediate between the parent species; 2) inheriting the 
trait of one of the parents 3) lacking the trait of one of the parents, or 4) 
showing a ‘new’ character absent from both parents. In the present study, the 
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TaBLE 1. Variables of Anemia organensis, Anemia villosa and the hybrid (A. organensis x A. 
villosa). Means+SD (standard deviation) and F values (ANOVA, *P<0,05). Means followed by 
different letters differ at P<0.05 (Tukey’s HSD test) 


A.organensis A. villosa hybrid F 
Epidermal thickness - adaxial 20.02+2.89a 23,06+5.55b 20.5143.52a 1457 1 
surface (jm) 
Epidermal thickness - abaxial 21.77+3.74a 25.82+7.09b 25.6343.33b 21.70* 
surface (1m) 
Stomatal length (um) 42.61+4.26a 52.67+4.34b 41.49+2.93a 191.81* 
Stomatal density (mm) 42.80+6.47a 47.83+3.56b 47.42+2.97b 123077 
Leaf thickness in vein region (j1m) 314.83+44.47a 219.05+29.22b 326.64+59.49a 132.18* 
Leaf thickness non-vein region (um) 506.72+51.38a 308.08+25.44b 357.52+42.33c 557.38* 
Mesophyll thickness (1m) 279.79+44.91a 16.57+15.73b 287.36457.37a  216.55* 
Trichome density - adaxial Ua garg laa EE 7.23+1.68b 321,064" -9192.425 
surface (mm) 
Trichome density - abaxial 15.17+3.45a 15.59+3.45a 26.90+4.35b 115.66* 
surface (mm) 
Number of trichome cells S'14+4,19a 6.80+0.94b 5.60+1.40c 3017" 
Spore length (m) 92.02+5.49a 84.18+6.23b 81.73+10.20b 55.39* 


hybrid has a vein on the leaf margin consisting only of xylem cells not present 
in the parent plants. With respect to quantitative variables, the hybrid differed 
from the parent plants in 2 of the 11 variables, with 5 variables similar to A. 
organensis and 4 to A. villosa. Thus, in relation to quantitative variables, the 
hybrid exhibits traits from both parent plants. 

Webb and Carlquist (1964) investigated Salvia species and found that 
trichome density on the abaxial and adaxial leaf surfaces was useful in 
discriminating hybrids from parents. In the present study, trichome density 
was higher on the adaxial surface in A. villosa and on the abaxial surface in the 
hybrid. 

There were significant differences in stomatal length between the taxa 
studied. The stomata of A. organensis and the hybrid showed similar values, 
but were larger in A. villosa. Similar data were obtained by Brandham and 
Cutler (1978) in a study with Aloe species. 

Multivariate exploratory techniques (PCA and DFA) were important in 
distinguishing between the hybrid and parent plants. The results showed 


TABLE 2. Factor loadings in PCA and factor structure coefficients in DFA for the variables of 
Anemia organensis, Anemia villosa and the hybrid (A. organensis x A. villosa). 


Variables PC1 PG? DF1 DF2 
Stomatal density 0.272 0.393 -0.406 0.526 
Thickness -0.903 0.120 0.917 -0.072 
Mesophyll thickness -0.907 0.101 0.055 -0.490 
Epidermal thickness - adaxial surface -0.219 0.164 -0.109 -0.038 
Trichome density - adaxial surface 0.729 -0.202 -0.021 0.468 
Trichome density - abaxial surface -0.361 0.785 -0.316 -0.432 
Spore length -0.205 -0.693 0.191 -0.084 
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Fic. 4. Discriminant analysis of Anemia villosa, Anemia organensis, and hybrid. 


100% separation between the groups. Thébaud and Abbot (1995), Gurgeli 
(1997), Levin (1999), and Albaladejo and Silvestre (2004) used the same 
techniques to separate hybrids, and the last authors also obtained 100% 
separation. 

In conclusion, hybrids can be distinguished from parent plants by their 
morphological and anatomical characteristics. Hybrid morphology is interme- 
diate in relation to parent plants, whereas the anatomy of hybrids shares 
features with parent plants and exhibits unique characteristics (leaf thickness 
and number of trichome cells). The hybrid is similar to A. organensis in five 
variables (adaxial surface epidermal thickness, stomatal length, leaf thickness 
in the non-vein region, mesophyll thickness, and trichome density on the 
adaxial surface) and to A. villosa in four variables (abaxial surface epidermal 
thickness, stomatal density, trichome density on the abaxial surface, and spore 


length). 
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ABstRACT.—We present here updates to the known species of ferns and lycophytes in the Hawaiian 
Islands since the publication of Palmer (2003) and the overview and updates of Vernon and Ranker 
(2013). 


Palmer (2003) represented the first thorough and concise presentation of the 
diversity of ferns and lycophytes on the Hawaiian Islands. Vernon and Ranker 
(2013) provided updates to the status of known species up to that time, 
including nomenclatural changes and the addition of newly named species 
and new records of naturalized species. Since that last update one additional 
new species was described (Wood and Wagner, 2017) and several additional 
naturalized, non-native species have been documented (see references below). 
Thus, we present here an updated checklist of the native and naturalized ferns 
and lycophytes of the Hawaiian Islands based on the references cited below. 
The checklist follows the familial and generic classification system recom- 
mended by PPG I (2016). The checklist also includes known and putative 
hybrids and recognized forms, not included in Vernon and Ranker (2013). 

A spreadsheet version of this checklist may be downloaded at: http://www. 
botany. hawaii.edu/people/tom-ranker/ 

Abbreviations: = endemic; ! indigenous; naturalized; (ex)=extinct; 
USFWS=U. S. Fish & Wildlife Service; N=Ni‘ihau; Le=Lehua; K=Kaua‘i; 
O=O*‘ahu; Mo=Moloka‘i; L=Lana‘i; Ma=Maui; Ka=Kaho‘olawe; H=Hawai‘i 

* Synonyms listed are those that represent updates subsequent to Palmer 
(2003) j 


NZ 


LYCOPHYTES 


1. Lycopodiaceae P.Beauv. ex Mirb. 
Huperzia Bernh. 

" Huperzia erosa Beitel & W.H.Wagner 
Distribution: K/O/Mo/L/Ma/H 

Huperzia erubescens (Brack.) Holub 
Distribution: K/O/Mo/Ma/H 

: Huperzia haleakalae (Brack.) Holub 
Distribution: Ma(ex) 

Huperzia serrata (Thunb. ex Murray) Trevis. 
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Distribution: K/O/Mo/L/H 
" Huperzia subintegra (Hillebr.) Beitel & W.H.Wagner 
Distribution: K/O/Mo/Ma 
© Huperzia Xcarlquistii Beitel & W.H.Wagner (Huperzia erubescens X 
H. subintegra) 
Distribution: O/Ma 
5, Huperzia erubescens (Brack.) Holub X H. haleakalae (Brack.) Holub 
Distribution: K/Ma 
" Huperzia xgillettii Beitel & W.H.Wagner (Huperzia serrata < H. 
subintegra) 
Distribution: K/O/Mo/Ma/H 
© Huperzia <medeirosii Beitel & W.H.Wagner (Huperzia haleakalae 
H. subintegra) 
Distribution: Ma 
" Huperzia sulcinervia (Spring) Trevis. (Huperzia erubescens < H. 
serrata) 
Distribution: O/L/H 
Lycopodium L. 
' Lycopodium venustulum Gaudich. var. venustulum 
Distribution: K/O/Mo/L/Ma/H 
© Lycopodium venustulum Gaudich. var. verticale W.H.Wagner 
Distribution: H 
Palhinhaea Franco & Carv. 
' Palhinhaea cernua (L.) Franco & Vasc. 
Syn. Lycopodiella cernua (L.) Pic.Serm. 
Distribution: K/O/Mo/L/Ma/H 
Phlegmariurus (Herter) Holub 
' Phlegmariurus filiformis (Sw.) W.H.Wagner 
Syn. Huperzia filiformis (Sw.) Holub 
Distribution: K/O/Mo/L/Ma/H 
© Phlegmariurus mannii (Hillebr.) W.H.Wagner 
Syn. Huperzia mannii (Hillebr.) Kartesz & Gandhi 
Distribution: K(ex)/Ma/H 
USFWS category: Endangered 
: Phlegmariurus nutans (Brack.) W.H.Wagner 
Syn. Huperzia nutans (Brack.) Rothm. 
Distribution: K(ex)/O 
USFWS category: Endangered 
' Phlegmariurus phyllanthus (Hook. & Arn.) R.D.Dixit 
Syn. Huperzia phyllantha (Hook. & Arn.) Holub 
Distribution: K/O/Mo/L/Ma/H 
© Phlegmariurus stemmermanniae A.C.Medeiros & W.H.Wagner 
Syn. Huperzia stemmermanniae (A.C.Medeiros & W.H.Wagner) 
Kartesz 
Distribution: Ma/H 
USFWS category: Endangered 
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© Phlegmariurus <koolauensis W.H.Wagner (Phlegmariurus nutans x 

P. phyllanthus) 
Distribution: O 

2. Isoétaceae Reichenb. 
Isoétes L. 

© Isoétes hawaiiensis W.C.Taylor & W.H.Wagner 

Distribution: Ma/H 
3. Selaginellaceae Willk. 
Selaginella P. Beauv. 

' Selaginella arbuscula (Kaulf.) Spring 
Distribution: K/O/Mo/L/Ma/H 
Note: Selaginella arbuscula was previously considered an endemic 
species. This species has a greater distribution, occurring in the 
Society Islands, Ualan, Santa Cruz Island (Vanikoro), and the 
Marquesas Islands (Hassler and Swale, 2019). 

® Selaginella deflexa Brack. 
Distribution: K/O/Mo/Ma/H 

NZ Selaginella kraussiana (Kunze) A.Braun 
Distribution: O/Ma/H 

N% Selaginella laxa Spring 
Distribution: K 
New record based on specimens: Palmer 582 & Palmer 583 (HAW) 

N% Selaginella stellata Spring 
Distribution: H 

N% Selaginella umbrosa Lem. ex Hieron. 
Distribution: H 


FERNS 


1. Ophioglossaceae Martinov 
Ophioderma (Blume) Endl. 
' Ophioderma pendulum (L.) G.Pres| 
subsp. falcatum (C.Presl) R.T.Clausen 
Distribution: K/O/Mo/L/Ma/H 
subsp. pendulum 
Distribution: O/Mo/Ma/H 
Ophioglossum L. 
' Ophioglossum nudicaule L.f. 
Distribution: K/O/H 
' Ophioglossum petiolatum Hook. 
_ Distribution: K/O/Mo/L/Ma/H 
‘ Ophioglossum polyphyllum A.Braun 
Distribution: K/O/Mo/L(ex)/Ka/Ma/H 
Sceptridium Lyon 
" Sceptridium subbifoliatum (Brack.) Lyon 
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Distribution: K(ex)/O(ex)/Mo(ex)/L(ex)/Ma(ex)/H(ex) 
2. Psilotaceae J.W.Griff. & Henfr. 
Psilotum Sw. 
' Psilotum complanatum Sw. 
Distribution: K/O/Mo/L/Ma/H 
' Psilotum nudum (L.) P. Beauv. 
Distribution: N/K/O/Mo/L/Ka/Ma/H 
' Psilotum Xintermedium W.H.Wagner (Psilotum complanatum  P. 
nudum) 
Distribution: K/O/L/Ma 
3. Marattiaceae Kaulf. 
Angiopteris Hoffm. 
N% Angiopteris evecta (G.Forst.) Hoffm. 
Distribution: K/O/Mo/L/Ma/H 
Marattia Sw. 
E Marattia douglasii (C.Presl) Baker 
Distribution: K/O/Mo/L/Ma/H 
4. Hymenophyllaceae Mart. 
Callistopteris Copel. 
" Callistopteris baldwinii (D.C.Eaton) Copel. 
Distribution: K/O/Mo/L/Ma/H 
Crepidomanes C. Presl 
© Crepidomanes draytonianum (Brack.) Ebihara & K.Iwats. 
Syn. Vandenboschia draytoniana (Brack.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
: Crepidomanes minutum (Blume) K.Iwats. 
Syn. Gonocormus minutus (Blume) Bosch; Crepidomanes proliferum 
(Blume) Bostock; Gonocormus prolifer (Blume) Prantl 
Distribution: K/O/Mo/L/Ma/H 
Hymenophyllum Sm. 
© Hymenophyllum Ianceolatum Hook. & Arn. 
Syn. Sphaerocionium lanceolatum (Hook. & Arn.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
© Hymenophyllum obtusum Hook. & Arn. 
Syn. Sphaerocionium obtusum (Hook. & Arn.) Copel. 
Distribution: O/Mo/L/Ma/H 
© Hymenophyllum recurvum Gaudich. 
Syn. Mecodium recurvum (Gaudich.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
Vandenboschia Copel. 
© Vandenboschia cyrtotheca (Hillebr.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
E Vandenboschia davallioides (Gaudich.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
© Vandenboschia tubiflora F.S.Wagner 
Distribution: K 
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5. Gleicheniaceae C. Pres] 
Dicranopteris Bernh. 
Dicranopteris linearis (Burm.f.) Underw. 
Distribution: K/O/Mo/L/Ma/H 
Includes: D. linearis (Burm. f.) Underw. f. linearis and D. linearis 
(Burm. f.) Underw. f. emarginata (Brack.) W.H.Wagner 
Diplopterygium (Diels) Nakai 
" Diplopterygium pinnatum (Kunze) Nakai 
Distribution: K/O/Mo/L/Ma/H 
Sticherus C. Pres] 
© Sticherus owhyhensis (Hook.) Ching 
Distribution: K/O/Mo/L/Ma/H 
6. Lygodiaceae M.Roem. 
Lygodium Sw. 
N“ Lygodium japonicum (Thunb.) Sw. 
Distribution: O/H 
7. Schizaeaceae Kaulf. 
Schizaea Sm. 
E Schizaea robusta Baker 
Distribution: K/O/Mo/L/Ma/H 
8. Marsileaceae Mirb. 
Marsilea L. 
NZ Marsilea minuta L. 
Syn. Marsilea crenata C.Presl 
Distribution: O 
E Marsilea villosa Kaulf. 
Distribution: N/O/Mo 
USFWS category: Endangered 
9. Salviniaceae Martinov 
Azolla Lam. 
N¢ Azolla filiculoides Lam. 
Distribution: K/O/Mo/L/Ma/H 
Salvinia Seg. 
NZ Salvinia molesta D.S.Mitch. 
Distribution: O/H 
10. Cibotiaceae Korall 
Cibotium Kaulf. 
= Cibotium chamissoi Kaulf. 
Distribution: O/Mo/L/Ma/H 
" Cibotium glaucum (Sm.) Hook. & Arn. 
Distribution: K/O/Mo/L/Ma/H 
" Cibotium menziesii Hook. 
Distribution: K/O/Mo/L/Ma/H 
" Cibotium nealiae O.Deg. 
Distribution: K 
" Cibotium xheleniae D.D.Palmer (Cibotium chamissoi X C. menziesii) 
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Distribution: O 
11. Cyatheaceae Kaulf. 
Sphaeropteris Bernh. 
N Sphaeropteris cooperi (Hook. ex F.Muell.) R.M.Tryon 
Distribution: K/O/L/Ma/H 
12. Dicksoniaceae M.R. Schomb. 
Dicksonia L’Her. 
N% Dicksonia fibrosa Colenso 
Distribution: H 
New record (Lorence and Flynn, 2006) 
13. Lindsaeaceae C. Pres] ex M.R. Schomb. 
Lindsaea Dryand. ex Sm. 
N¢ Lindsaea ensifolia Sw. 
Distribution: K/O/Mo/Ma/H 
Lindsaea repens (Bory) Thwaites 
E var. macraeana (Hook. & Arn.) G.Chr. 
Distribution: K/O/Mo/L/Ma/H 
Odontosoria Fée 
' Odontosoria chinensis (L.) J.Sm. 
Syn. Sphenomeris chinensis (L.) Maxon 
Distribution: K/O/Mo/L/Ma/H 
xLindsaeosoria W.H.Wagner 
© xLindsaeosoria flynnii W.H.Wagner (Lindsaea ensifolia X Odonto- 
soria chinensis) 
Distribution: K/O 
New island record for O (Frohlich and Lau, 2014) 
14. Dennstaedtiaceae Lotsy 
Hypolepis Bernh. 
© Hypolepis hawaiiensis Brownsey 
var. hawaiiensis 
Distribution: K/O/Mo/L/Ma/H 
var. mauiensis (Hillebr.) D.D.Palmer 
Distribution: Ma 
USFWS category: Endangered 
Microlepia C. Pres] 
' Microlepia speluncae (L.) T.Moore 
Distribution: K/O/Mo/Ma/H 
Microlepia strigosa (Thunb.) C.Presl 
E var. mauiensis (W.H.Wagner) D.D.Palmer 
Distribution: O/Ma/H 
USFWS category: Endangered 
' var. strigosa 
Distribution: K/O/Mo/L/Ma/H 
Microlepia Xadulterina W.H.Wagner (Microlepia speluncae x M. 
strigosa) 
Distribution: K/O/H 
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Pteridium Gled. ex Scop. 
Pteridium aquilinum (L.) Kuhn 
S subsp. decompositum (Gaudich.) Lamoureux ex J.A. Thomson 
Distribution: K/O/Mo/L/Ma/H 
15. Pteridaceae E.D.M.Kirchn. 
Adiantum L. 
N% Adiantum ‘Edwinii’ 
Distribution: O/L/Ma 
Palmer (2003) noted that this taxon may be a cultivar of either A. 
raddianum C.Presl or A. concinnum Willd. or a hybrid Daley them. 
' Adiantum capillus-veneris L. 
Distribution: N/K/O/Mo/L/Ma/H 
NZ Adiantum hispidulum Sw. 
Distribution: K/O/Mo/L/Ka/Ma/H 
N@ Adiantum macrophyllum Sw. 
Distribution: O 
New record (Lau and Frohlich, 2015) 
N@ Adiantum raddianum C.Pres| 
Distribution: K/O/Mo/L/Ma/H 
N@ Adiantum tenerum Sw. 
Distribution: O/Ma 
Ceratopteris Brongn. 
N¢ Ceratopteris gaudichaudii Brongn. var. vulgaris Masuyama & 
Watano 
Distribution: O 
N¢ Ceratopteris thalictroides (L.) Brongn. 
Distribution: K/O/H 
Based on specimen annotations at BISH by S. Masuyama, November 
2018 
Cheilanthes Sw. 
NZ Cheilanthes viridis (Forssk.) Sw 
Distribution: K/O/L/Ma/H 
Coniogramme Fée 
© Coniogramme pilosa (Brack.) Hieron. 
Distribution: K/O/Mo/L/Ma/H 
Doryopteris J.Sm. 
" Doryopteris angelica K.R.Wood & W.H.Wagner 
Distribution: K 
USFWS category: Endangered 
" Doryopteris decipiens (Hook.) J.Sm. 
Distribution: N/K/O/Mo/L/Ka/Ma/H 
" Doryopteris decora Brack. 
Distribution: K/O/Mo/L/Ka/Ma/H 
" Doryopteris subdecipiens W.H.Wagner (Doryopteris decipiens X D. 
decora) 
Distribution: K/O/L/M/Ka/H 
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" Doryopteris takeuchii (W.H.Wagner) W.H.Wagner 
Distribution: O 
USFWS category: Endangered 
© Doryopteris tryonii O.Deg. & I.Deg. 
Distribution: K 
Note: Possibly a variant of D. decora; needs further study. 
See: Yesilyurt (2005) 
Haplopteris C. Pres] 
' Haplopteris elongata (Sw.) E.H.Crane 
Distribution: K/O/Mo/L/Ma/H 
Pellaea Link 
' Pellaea ternifolia (Cav.) Link 
Distribution: K/O/Mo/L/Ma/H 
Pityrogramma Link 
N¢ Pityrogramma austroamericana Domin 
Distribution: N/K/O/Mo/L/Ka/Ma/H 
N¢ Pityrogramma calomelanos (L.) Link 
Distribution: N/K/O/Mo/L/Ka/Ma/H 
© Pityrogramma <mckenneyi W.H.Wagner (Pityrogramma austroamer- 
icana X P. calomelanos) 
Distribution: O 
Pteris L. 
' Pteris cretica L. 
Distribution: K/O/Mo/L/Ma/H 
© Pteris hillebrandii Copel. 
Distribution: K/O/Mo/L/Ma/H 
© Pteris irregularis Kaulf. 
Distribution: K/O/Mo/L/Ma/H 
E Pteris lidgatei (Baker) Christ 
Distribution: O/Mo/Ma 
USFWS category: Endangered 
' Pteris terminalis Wallich ex J.Agardh var. terminalis 
Syn. Pteris excelsa Gaudich. 
Distribution: K/O/Mo/L/Ma/H 
Note: See Ebihara et al. (2017) for explanation of synonym 
N¢ Pteris tremula R.Br. 
Distribution: Ma 
New record (Oppenheimer, 2007) 
N¢ Pteris vittata L. 
Distribution: K/O/Mo/L/Ma/H 
16. Cystopteridaceae (Payer) Shmakov 
Cystopteris Bernh. 
" Cystopteris douglasii Hook. 
Distribution: Ma/H 
© Cystopteris sandwicensis Brack. 
Distribution: K/O/L/Ma 
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17. Aspleniaceae Newman 
Asplenium L. 
© Asplenium acuminatum Hook. & Arn. 
Distribution: K/O/Mo/L/Ma/H 
' Asplenium adiantum-nigrum L. 
Distribution: K/O/Mo/L/Ma/H 
| Asplenium aethiopicum (Burm.f.) Bech. 
Distribution: K/O/Mo/L/Ma/H 
' Asplenium caudatum G.Forst. 
Syn. Asplenium horridum Kaulf. var. glabratum (W.J.Rob.) 
D.D.Palmer; Asplenium horridum Kaulf. var. horridum 
Distribution: K/O/Mo/L/Ma/H 
Asplenium contiguum Kaulf. 
' var. contiguum 
Distribution: K/O/Mo/L/Ma/H 
Note: Species listed as endemic in Palmer (2003), but reported also in 
Vanuatu by National Museum of Nature and Science (2008) and 
Ranker pers. obs. but needs critical evaluation. 
© var. hirtulum C.Chr. 
Distribution: K/Ma 
Asplenium dielerectum Viane 
Syn. Diellia erecta Brack. 
Distribution: O/Mo/L(ex)/Ma/H 
Asplenium dielfalcatum Viane 
Syn. Diellia falcata Brack. 
Distribution: O 
USFWS category: Endangered 
Asplenium diellaciniatum Viane 
Distribution: K 
USFWS category: Endangered 
See Lorence et al. (2013) for a discussion of this species. 
Asplenium dielmannii Viane 
Syn. Diellia mannii (D.C.Eaton) W.J.Rob. 
Distribution: K 
USFWS category: Endangered 
Asplenium dielpallidum N.Snow 
Syn. Diellia pallida W.H.Wagner 
Distribution: K 
USFWS category: Endangered 
Asplenium haleakalense W.H.Wagner 
Distribution: Ma/H 
Asplenium hobdyi W.H.Wagner 
Distribution: K/Mo/Ma/H 
' Asplenium insiticium Brack. __ 
_ Distribution: K/O/Mo/L/Ma/H 
" Asplenium kaulfussii Schltdl. 
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Includes: A. kaulfussii Schltdl. f. kaulfussii, A. kaulfussii Schltdl. f. 
gemmiparum (Hillebr.) D.D.Palmer, A. kaulfussii Schltdl. f. bipinna- 
tum (Hillebr.) D.D.Palmer, and A. kaulfussii Schltdl. f. dareoides 
(Hillebr.) W.H.Wagner 
Distribution: K/O/Mo/L/Ma/H 
© Asplenium leucostegioides Baker 
Syn. Diellia leucostegioides (Baker) W.H.Wagner 
Distribution: Ma(ex) 
' Asplenium lobulatum Mett. 
Distribution: K/O/Mo/L/Ma/H 
' Asplenium macraei Hook. & Grev. 
Status changed from Endemic to Indigenous; discovered in Marquesas 
Islands; see checklist: https://naturalhistory2.si.edu/botany/ 
marquesasflora/speciesdescr.cfm?genus=Asplenium&species= 
macraei; also specimens from Savii, Samoa: Lorence et al. 10572 
(PTBG, UC) and Flynn et al. 8626 (K, MBK, NY, PTBG, UC, US) 
Distribution: K/O/Mo/L/Ma/H 
' Asplenium monanthes L. 
Distribution: Ma/H 
' Asplenium nidus L. 
Distribution: K/O/Mo/L/Ma/H 
' Asplenium normale D.Don 
Distribution: K/O/Mo/L/Ma/H 
Asplenium peruvianum Desv. 
E var. insulare (C.V.Morton) D.D.Palmer 
Distribution: Ma/H 
USFWS category: Endangered 
' Asplenium polyodon G.Forst. 
Distribution: K/O/Mo/L/Ma/H 
© Asplenium schizophyllum C.Chr. 
Distribution: K/H 
© Asplenium sphenotomum Hillebr. 
Distribution: K/O/Mo/L/Ma/H 
Asplenium trichomanes L. 
© subsp. densum (Brack.) W.H.Wagner 
Distribution: Ma/H 
ss Asplenium unisorum (W.H.Wagner) Viane 
Syn. Diellia unisora W.H.Wagner 
Distribution: O 
USFWS category: Endangered 
 Asplenium <flagrum W.H.Wagner & D.D.Palmer (Asplenium hobdyi x 
A. normale) 
Distribution: K/M 
New island record for M (Oppenheimer, 2008) 
© Asplenium joellauii N.Snow (Asplenium dielfalcatum X A. unisorum) 
Distribution: O 
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New combination (Snow, 2011) 
© Asplenium <kokeense W.H.Wagner (Asplenium aethiopicum X A. 
polyodon) 
Distribution: K/O 
" Asplenium Xsphenocookii W.H.Wagner (Asplenium polyodon X A. 
sphenotomum) 
Distribution: K 
© Asplenium <waikamoi W.H.Wagner (Asplenium acuminatum X A. 
aethiopicum) 
Distribution: Ma 
Hymenasplenium Hayata 
} Hymenasplenium excisum (C.Presl) S.Linds. 
Syn. Asplenium excisum C.Pres] 
Distribution: K/O/Mo/L/Ma/H 
' Hymenasplenium unilaterale (Lam.) Hayata 
Syn. Asplenium unilaterale Lam. 
Distribution: K/O/Mo/L/Ma/H 
18. Thelypteridaceae Ching ex Pic.Serm. 
Amauropelta Kunze 
© Amauropelta globulifera (Brack.) Holttum 
Syn. Thelypteris globulifera (Brack.) C.F.Reed 
Distribution: K/O/Mo/L/Ma/H 
Christella H.Lév. 
© Christella boydiae (D.C. Eaton) Holttum 
Syn. Cyclosorus boydiae (D. C. Eaton) W.H.Wagner 
Distribution: O/Ma 
= Christella cyatheoides (Kaulf.) Holttum 
Syn. Cyclosorus cyatheoides (Kaulf.) Farw. 
Distribution: K/O/Mo/L/Ma/H 
N’ Christella dentata (Forssk.) Brownsey & Jermy 
Syn. Cyclosorus dentatus (Forssk.) Ching 
Distribution: N/K/O/Mo/L/Ma/H 
NZ Christella parasitica (L.) H.Lév. 
Syn. Cyclosorus parasiticus (L.) Farw. 
Distribution: K/O/Mo/L/Ma/H 
" Christella wailele (Flynn) D.D.Palmer 
Syn. Cyclosorus wailele (Flynn) W.H.Wagner 
Distribution: K 
§  Christella cyatheoides (Kaulf.) Holttum Xx C. dentata (Forssk.) 
Brownsey & Jermy 
Syn. Cyclosorus Xpalmeri (W.H.Wagner) W.H.Wagner; Thelypteris 
<palmeri W.H.Wagner 
Distribution: K/O . 
" Christella dentata (Forssk.) Brownsey & Jermy X C. parasitica (L.) 
H.Lév. 
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Syn. Christella Xincesta (W.H.Wagner) Nakaike & Kawabata; Chris- 
tella intermedius (W.C.Shieh & J.L.Tsai) D.D.Palmer 
Distribution: K/O/Ma 
Cyclosorus Link 
, Cyclosorus interruptus (Willd.) H.It6 
Distribution: K/O/Mo/L/Ma/H 
Macrothelypteris (H. It6) Ching 
N¢ Macrothelypteris torresiana (Gaudich.) Ching 
Distribution: K/O/Mo/Ma/H 
New island record for Mo (Oppenheimer, 2016) 
Pneumatopteris Nakai 
= Pneumatopteris hudsoniana (Brack.) Holttum 
Syn. Cyclosorus hudsonianus (Brack.) Ching 
Distribution: K/O/Mo/L/Ma/H 
© Pneumatopteris pendens D.D.Palmer 
Syn. Cyclosorus pendens (D.D.Palmer) N.Snow 
Distribution: K/O/Mo/Ma/H 
x Pneumatopteris sandwicensis (Brack.) Holttum 
Syn. Cyclosorus sandwicensis (Brack.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
Pseudophegopteris Ching 
© Pseudophegopteris keraudreniana (Gaudich.) Holttum 
Distribution: K/O/Mo/L/Ma/H 
19. Blechnaceae Newman 
Blechnopsis C.Pres] 
Ne Blechnopsis orientalis (L.) C.Pres] 
Syn. Blechnum orientale L. 
Distribution: O 
New State record (Lau and Frohlich, 2012) 
Blechnum L. 
N% Blechnum appendiculatum Willd. 
Distribution: K/O/Mo/L/Ma/H 
Doodia R. Br. 
E Doodia kunthiana Gaudich. 
Distribution: K/O/Mo/L/Ma/H 
© Doodia lyonii O. Deg. 
Distribution: K/O/Ma/H(ex) 
Sadleria Kaulf. 
© Sadleria cyatheoides Kaulf. 
Distribution: K/O/Mo/L/Ma/H 
© Sadleria pallida Hook. & Arn. 
Distribution: K/O/Mo/L/Ma/H 
E Sadleria souleyetiana (Gaudich.) T.Moore 
Distribution: K/O/Mo/L/Ma/H 
© Sadleria squarrosa (Gaudich.) T.Moore 
Distribution: K/O/Mo/L/Ma/H 
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© Sadleria unisora (Baker) W.J.Rob. 
Distribution: K 
© Sadleria wagneriana D.D.Palmer & Flynn 
Distribution: K 
20. Athyriaceae Alston 
Athyrium Roth 
© Athyrium haleakalae K.R.Wood & W.L.Wagner 
Distribution: Ma 
Newly described species (Wood and Wagner, 2017) 
© Athyrium microphyllum (Sm.) Alston 
Distribution: K/O/Mo/L/Ma/H 
Deparia Hook. & Grev. 
© Deparia cataracticola M.Kato 
Distribution: K 
: Deparia fenzliana (Luerss.) M.Kato 
Distribution: K/O/Mo/L/Ma/H 
© Deparia kaalaana (Copel.) M.Kato 
Distribution: K(ex)/Ma(ex)/H(ex) 
USFWS category: Endangered 
© Deparia marginalis (Hillebr.) M.Kato 
Distribution: K/Mo/L/Ma/H 
N¢ Deparia petersenii (Kunze) M.Kato subsp. petersenii 
Distribution: K/O/Mo/L/Ma/H 
© Deparia prolifera (Kaulf.) Hook. & Grev. 
Distribution: K/O/Mo/L/Ma/H 
Diplazium Sw. 
© Diplazium arnottii Brack. 
Distribution: K/O/Mo/L/Ma/H 
re Diplazium esculentum (Retz.) Sw. 
Distribution: K/O/Mo/L/Ma/H 
" Diplazium molokaiense W.J.Rob. 
Distribution: K(ex)/O(ex)/Mo(ex)/L(ex)/Ma 
USFWS category: Endangered 
a Diplazium sandwichianum (C. Presl) Diels 
Distribution: K/O/Mo/L/Ma/H 
21. Dryopteridaceae Herter 
Arachniodes Blume 
© Arachniodes insularis W.H.Wagner 
Distribution: K/O/Mo/Ma/H 
Ctenitis (C. Chr.) C. Chr. 
" Ctenitis latifrons (Brack.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
© Ctenitis squamigera (Hook. & Arn.) Copel. 
Distribution: K/O/Mo/L/Ma 
USFWS category: Endangered 
Cyrtomium C. Presl 
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, Cyrtomium caryotideum (Wall. ex Hook. & Grev.) C.Pres| 
Distribution: K/O/Mo/L/Ma/H 
pe Cyrtomium falcatum (L.f.) C.Pres] 
Distribution: K/O/Mo/L/Ma/H 
Dryopteris Adans. 
s Dryopteris crinalis (Hook. & Arn.) C.Chr. 
var. crinalis 
Distribution: K/O/Mo/L/Ma/H 
var. podosora (W.H.Wagner) D.D.Palmer 
Distribution: K 
USFWS category: Endangered 
© Dryopteris fuscoatra (Hillebr.) W.J.Rob. 
var. fuscoatra 
Distribution: K/O/Mo/L/Ma/H 
var. Jamoureuxii Fraser-Jenk. 
Distribution: Ma 
Dryopteris glabra (Brack.) Kuntze 
var. alboviridis (W.H.Wagner) D.D.Palmer 
Distribution: K 
var. flynnii D.D.Palmer 
Distribution: K 
var. glabra 
Distribution: K/O/Mo/L/Ma/H 
var. hobdyana (W.H.Wagner) D.D.Palmer 
Distribution: Ma 
var. nuda (Underw.) Fraser-Jenk. 
Distribution: K/O/Mo/Ma 
var. pusilla (Hillebr.) Fraser-Jenk. 
Distribution: K 
USFWS category: Endangered 
var. soripes (Hillebr.) Herat ex Fraser-Jenk. 
Distribution: K/O/Mo/Ma/H 
© Dryopteris hawaiiensis (Hillebr.) W.J.Rob. 
Distribution: K/O/Mo/Ma/H 
© Dryopteris mauiensis C. Chr. 
Distribution: K/O/Mo/L/Ma/H 
© Dryopteris rubiginosa (Brack.) Kuntze 
Syn. Nothoperanema rubiginosum (Brack.) A.R.Sm. & D.D.Palmer 
Distribution: K/O/Mo/L/Ma/H 
" Dryopteris sandwicensis (Hook. & Arn.) C.Chr. 
Distribution: K/O/Mo/L/Ma/H 
© Dryopteris subbipinnata W.H.Wagner & Hobdy 
Distribution: Ma/H 
© Dryopteris tetrapinnata W.H.Wagner & Hobdy 
Distribution: Ma 
© Dryopteris unidentata (Hook. & Arn.) C.Chr. 
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var. paleacea (Hillebr.) Herat ex Fraser-Jenk. 
Distribution: K/O/Mo/Ma/H 
var. unidentata 
Distribution: K/O/Mo/L/Ma/H 
' Dryopteris wallichiana (Spreng.) Hy]. 
Distribution: K/O/Mo/Ma/H 
Elaphoglossum Schott ex J.Sm. 
© Elaphoglossum aemulum (Kaulf.) Brack. 
Distribution: K/O/Mo/L/Ma/H(ex) 
© Elaphoglossum alatum Gaudich. 
Distribution: O 
© Elaphoglossum crassicaule Copel. 
Distribution: K 
" Elaphoglossum crassifolium (Gaudich.) W.R.Anderson & Crosby 
Distribution: K/O/Mo/L/Ma/H 
© Elaphoglossum fauriei Copel. 
Distribution: O/Mo 
' Elaphoglossum paleaceum (Hook. & Grev.) Sledge 
Distribution: K/O/Mo/L/Ma/H 
© Elaphoglossum parvisquameum Skottsb. 
Distribution: Mo/L/Ma/H 
" Elaphoglossum pellucidum Gaudich. 
Distribution: K/O/Mo/Ma/H 
© Elaphoglossum wawrae (Luerss.) C.Chr. 
Distribution: K/O/Mo/Ma/H 
Elaphoglossum alatum Gaudich. x E. fauriei Copel. 
Distribution: O 
" Elaphglossum fauriei Copel. x E. parvisquameum Skottsb. 
Distribution: Mo 
Polystichum Roth 
" Polystichum bonseyi W.H.Wagner & Hobdy 
Distribution: Ma/H 
" Polystichum haleakalense Brack. 
Distribution: Ma/H 
" Polystichum hillebrandii Carruth. 
Distribution: Ma/H 
22. Nephrolepidaceae Pic.Serm. 
Nephrolepis Schott 
N% Nephrolepis brownii (Desv.) Hovenkamp & Miyam. 
Syn. Nephrolepis multiflora (Roxb.) F.M.Jarrett ex C.V.Morton 
Distribution: N/K/O/Mo/L/Ka/Ma/H 
' Nephrolepis cordifolia (L.) C.Presl 
Distribution: K/O/Mo/L/Ma/H 
Nephrolepis exaltata (L.) Schott — 
" subsp. hawaiiensis W.H.Wagner 
Distribution: N/K/O/Mo/L/Ma/H 
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N“ Nephrolepis falcata (Cav.) C.Chr. ‘Furcans’ 
Distribution: K/O/Mo/L/Ma/H 
N* Nephrolepis hirsutula (G. Forst.) C.Presl ‘Superba’ 
Distribution: K/Ma 
" Nephrolepis <copelandii W.H.Wagner (Nephrolepis cordifolia x N. 
brownii) 
Distribution: O 
© Nephrolepis Xmedlerae W.H.Wagner (Nephrolepis exaltata subsp. 
hawaiiensis X N. brownii) 
Distribution: O 
23. Tectariaceae Panigrahi 
Tectaria Cav. 
E Tectaria gaudichaudii (Mett.) Maxon 
Distribution: K/O/Mo/L/Ma/H 
N¢ Tectaria incisa Cav. 
Distribution: K/O/Ma/H 
24. Davalliaceae M.R.Schomb. ex A.B.Frank 
Davallia Sm. 
N% Davallia fejeensis Hook. 
Distribution: O 
25. Polypodiaceae Bercht. & J.Pres| 
Adenophorus Gaudich. 
5 Adenophorus abietinus (D. C. Eaton) K.A.Wilson 
Distribution: K/O/L/Ma 
Adenophorus epigaeus (L.E.Bishop) W.H.Wagner 
Distribution: K 
Adenophorus haalilioanus (Brack.) K.A.Wilson 
Distribution: K/O 
Adenophorus hymenophylloides (Kaulf.) Hook. & Grev. 
Distribution: K/O/Mo/L/Ma/H 
Adenophorus oahuensis (Copel.) L.E.Bishop 
Distribution: O 
Adenophorus periens L.E.Bishop 
Distribution: K/O(ex)/Mo/L(ex)/Ma(ex)/H 
USFWS category: Endangered 
Note: Extant population status on K, Mo, and H needs verification - 
possibly extinct 
Adenophorus pinnatifidus Gaudich. 
var. pinnatifidus 
Distribution: K/O/Mo/L/Ma/H 
var. rockii (Copel.) D.D.Palmer 
Distribution: K/O/Mo/Ma 
Adenophorus tamariscinus (Kaulf.) Hook. & Grev. 
var. montanus (Hillebr.) L.E.Bishop 
Distribution: Mo/Ma/H 
var. tamariscinus 
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Distribution: K/O/Mo/L/Ma/H 
" Adenophorus tenellus (Kaulf.) Ranker 
Syn. Gramumitis tenella Kaulf. 
Distribution: K/O/Mo/L/Ma/H 
© Adenophorus tripinnatifidus Gaudich. 
Distribution: K/O/Mo/L/Ma/H 
© Adenophorus Xabbottiae W.H.Wagner (Adenophorus hymenophyl- 
loides X A. tamariscinus) 
Distribution: K/O 
© Adenophorus Xcarsonii Ranker (Adenophorus hy m6 Oe x 
A. tripinnatifidus) 
Distribution: H 
" Adenophorus oahuensis (Copel.) L.E.Bishop < A. pinnatifidus 
Gaudich. 
Syn. Oligadenus Xbishopii W.H.Wagner 
Distribution: O 
Aglaomorpha Schott. 
N@ Aglaomorpha rigidula (Sw.) Hovenkamp & S.Linds. 
Syn. Drynaria rigidula (Sw.) Bedd. 
Distribution: O 
New record (Lau and Frohlich, 2013; reported as D. rigidula) 
Lepisorus (J.Sm.) Ching 
' Lepisorus thunbergianus (Kaulf.) Ching 
Distribution: K/O/Mo/L/Ma/H 
Microsorum Link 
N@ Microsorum grossum (Langsd. & Fisch.) S.B.Andrews — 
Syn. Phymatosorus grossus (Langsd. & Fisch.) Brownlie 
Distribution: K/O/Mo/L/Ma/H 
Note: M. grossum and M. scolopendria may be conspecific; see 
Nooteboom (1997) and Krier et al. (2008) 
N* Microsorum scolopendria (Burm.f.) Copel. 
Syn. Phymatosorus scolopendria (Burm.f.) Pic.Serm. 
Distribution: Ma 
M. grossum and M. scolopendria may be conspecific; see Nooteboom 
(1997) and Krier et al. (2008) 
© Microsorum spectrum (Kaulf.) Copel. 
var. pentadactylum (Hillebr.) D.D.Palmer 
Distribution: K/Ma 
var. spectrum 
Distribution: K/O/Mo/L/Ma/H 
Oreogrammitis Copel. 
e Oreogrammitis baldwinii (Baker) Parris 
Syn. Grammitis baldwinii (Baker) Copel. 
Distribution: K ) 
" Oreogrammitis forbesiana (W.H.Wagner) Parris 
Syn. Grammitis forbesiana W.H.Wagner 
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Distribution: K/Mo/Ma 
s Oreogrammitis hookeri (Brack.) Parris 
Syn. Grammitis hookeri (Brack.) Copel. 
Distribution: K/O/Mo/L/Ma/H 
Phlebodium (R.Br.) J.Sm. 
NZ Phlebodium aureum (L.) J.Sm. 
Distribution: K/O/Mo/L/Ma/H 
Platycerium Desv. 
N¢ Platycerium bifurcatum (Cav.) C.Chr. 
Distribution: O/Ma/H 
N¢ Platycerium superbum de Jonch. & Hennipman 
Distribution: O 
Polypodium L. 
© Polypodium pellucidum Kaul. 
var. acuminatum D.D.Palmer 
Distribution: K/Ma 
var. pellucidum 
Distribution: K/O/Mo/L/Ma/H 
Includes: P. pellucidum Kaulf. var. pellucidum f. pellucidum and P. 
pellucidum Kaulf. var. pellucidum f. opacum (Hillebr.) D.D.Palmer 
var. vulcanicum Skottsb. 
Distribution: Mo/Ma/H 
Pyrrosia Mirb. 
N¢ Pyrrosia longifolia (Burm.f.) C.V.Morton 
Distribution: O 
New record (Frohlich and Lau, 2014) 
N¢ Pyrrosia piloselloides (L.) M.G.Price 
Distribution: O 
New record (Lau and Frohlich, 2012) 
Stenogrammitis Labiak 
5 Stenogrammitis saffordii (Maxon) Labiak 
Syn. Lellingeria saffordii (Maxon) A.R.Sm. & R.C.Moran 
Distribution: K/O/Mo/L/Ma/H 
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SHORTER NOTE 


Occasional Bifid Leaves in Selaginella helvetica.—The lycophyte leaf is 
unique having evolved independently of the leaves of euphyllophytes 
(Tomescu, Trends in Plant Science 14:5-12. 2008). The leaves of lycophytes 
are commonly referred to as microphylls and first formed as enations of the 
stem. Into this enation, a single vascular trace extending from the stele 
infiltrated to facilitate hydration and the transport of photosynthates (Stewart 
and Rothwell, Paleobotany and the evolution of plants Second Edition, 
Cambridge University Press, Cambridge. 1993). Unlike the leaves of euphyllo- 
phytes, the unique developmental origin of lycophyte microphylls is believed 
to have placed considerable developmental constraint on the evolution of 
compound leaves and reticulate venation in this lineage (Harrison, Rezvani 
and Lugdale, Development 134:881-889. 2007). With only a few exceptions, 
microphylls in extant species have a single vascular trace in a simple, unlobed 
lamina. Reported variants to this single-vein vascular arrangement in 
microphylls include an occasional two or more veins in the simple leaves 
(or leaves with small lobes at the apex) of the extant species Selaginella 
willdenowii (Desv.) Baker (Grambast and Rosello, Comptes rendus de 
Académie des Sciences 261: 5183-5186. 1965) and S. martensii Spring 
(Webster, American Fern Journal 60:1—6. 1970). While the most unusual extant 
exception within the lycophyte clade is the branching veins in the simple 
leaves of the species S. schaffneri Hieron. (an example of this unique vascular 
arrangement is shown in Fig. 1A) and S. adunca A.Braun ex Hieron. (Wagner, 
Beitel, Wagner, Science 218:793—794. 1982). A new and undescribed variant to 
the simple and single-veined microphyll leaves of lycophytes can be found on 
occasion in S. helvetica (L.) Spring (Fig. 1C and 1D). During work involving the 
close observation of stems and leaves of an individual of this species, collected 
in Southern Germany and subsequently cultivated in the greenhouses of 
Purdue University, Indiana, a number of bifid leaves were observed across 
numerous shoots. Two unique forms were observed. The first form found in 
the larger ventral leaves, had two equal sized lobes each containing a vein 
originating from a common vascular trace. These veins bifurcate at the base of 
the leaf and extend to the apex of each of the lobes (Fig. 1). The second form, 
observed in smaller dorsal leaves, had a short single vein extending to the 
terminus of the sinus separating two uneven sized lobes (Fig. 1). In both bifid 
leaf forms, a single and entire ligule was present (Fig. 1). Bifid leaves are not 
known in any extant lycophyte species; however, these variants in S. helvetica 
indicate that the genetic or developmental capacity for such an innovation 
exists in this lineage. Branched or lobed leaves in lycophytes have been 
described in a number of extinct genera and species including Estinnophyton, 
Protolepidodendron, Spencerites insignis and the genus Leclercqia (Stewart 
and Rothwell, 1993). Leclercgia species have the most comprehensively 
described and well-preserved leaves of these extinct, lobe-leaved lycophyte 
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Fic. 1. A. The simple leaf of Selaginella schaffneri displaying a unique branching and 
dichotomising venation arrangement (from herbarium specimen PH00023406). B. A common 
simple leaf of S. helvetica. C. A bifid ventral leaf of S. helvetica with veins extending to the apex of 
the leaf. D. A bifid dorsal leaf of S. helvetica with vein extending to the terminus of the sinus of the 
two lobes. All leaves were cleared with dilute sodium hypochlorite solution and stained with 
Toluidine Blue O, the leaves of S. helvetica were taken from the same plant, voucher accession 
PUL N18230. Bars = 500 um. | 


species; the leaves of Leclercqia species are divided at the tip into five, and up 
to seven, lobes (Bonamo, Banks and Grierson, Botanical Gazette 149:222 — 239. 
1988). In some Leclercgia species, lateral veins may have branched from a 
central main vascular trace to supply each of these lobes (Bonamo and 
Greirson, Botanical Society of America Miscellaneous Series 160:42. 1981). 
Other extinct leaf shape or venation variants in lycophytes include the two- 
veined leaves of Sigillaria (Stewart and Rothwell, 1993), and the sagittate, not 
lanceolate, leaves of Haskinsia (Bonamo, Banks and Grierson, 1988). These 
examples suggest that a complex leaf form and venation arrangement is 
developmentally possible in the lycophyte lineage, and has been for much of 
its history. The foliar and venation variants so far described in lycophytes 


SHORTER NOTES 75 


mirror (at a smaller scale) much of the anatomical variation of euphyllophyte 
leaves. This raises the important question of why, given the developmental 
capacity, has this variation not been acted upon by natural selection? This is 
especially perplexing considering the physiological benefits of increased 
photosynthetic capacity conferred by reticulate venation (Brodribb, Feild and 
Jordan, Plant Physiology 144:1890-1898. 2007). A potential explanation for 
this lack of adoption of euphyllophyte leaf modification broadly across the 
lycophyte lineage remains unknown.—Scotr A. M. McApam, Purdue Univer- 
sity, Department of Botany and Plant Pathology, Purdue Center for Plant 
Biology, 915 W. State St, West Lafayette, IN, 47907, USA, email: smcadam@ 
purdue.edu. 
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COVER CAPTION: The simple leaf of Selaginella schaffneri displaying a unique branching 


and dichotomising venation arrangement. Photo credit: Scott McAdam 


